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1 .  INTRODUCTION 


The  prediction  of  radar  and  optical  signatures  resulting  from 
atmospheric  nuclear  detonations  requires  reliable  fireball  entrainment 
and  mixing  models.  Essential  to  these  predictions  is  a  sound  knowledge  of 
the  temperature,  chemical  species,  and  velocity  fields  which  are  not 
provided  by  available  data  from  nuclear  tests.  The  objective  of  this 
program  was  to  investigate  turbulent  mixing  and  entrainment  of  air  in 
fireballs  by  means  of  a  subscale  laboratory  simulation  and  to  provide  a 
data  base  to  support  the  development  and  evaluation  of  theoretical  fireball 
models  and  codrs. 

Tie  development  of  fireballs  from  low  yield  (<100  kt)  bursts  at  low 
altitudes  (<30  km)  will  be  dominated  by  buoyancy  forces.  The  initial 
stage  of  these  fireballs  resulting  from  large  energy  release  within  the 
atmosphere  is  characterized  by  spherical  symmetry,  the  domination  of  radia¬ 
tive  energy  transport,  and  the  propagation  of  a  strong  spherical  shock  wave. 

On  the  order  of  a  second  f.fter  detonation,  the  fireball  comes  to  pressure 
equilibrium  with  the  ambient  air  and  because  of  its  low  density,  begins  to 
rise.  As  the  fireball  rises,  circulation  causes  a  toroid  to  be  formed. 

This  series  of  laboratory  experiments  which  have  been  performed  were 
designed  to  simulate  these  "buoyant  rise"  fireballs  after  pressure  equili¬ 
bration. 

This  report  provides  a  complete  summary  of  this  experimental  test 
program  which  was  conducted  during  the  period  from  August  1971  through 
September  1973.  Although  the  emphasis  in  this  report  is  on  the  experimental 
results  and  conclusions  drawn  from  those  results,  descriptions  of  the  test 
facility  and  discussions  of  the  data  acquisition,  interpretation  and  reduction 
procedures  are  also  included. 

Several  topics  were  investigated  during  the  course  of  this  experiment 
and  many  interesting  results  u~e  discussed  in  detail.  Mean  size/rise 
characteristics  and  particle  tracking  photographs  are  presented  to  illustrate 
gross  flow  features.  Laboratory  turbulent  results  were  compared  to 
dimensional  data  from  a  full  scale  nuclear  event  (through  pr /per  scaling) 
and  showed  '.'cry  gcod  agreement  Results  obtained  during  this  study 
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illustrated  a  marked  difference  between  turbulent  and  laminar  flow  regimes, 
and  thus  the  emphasis  in  this  report  is  placed  on  the  turbulent  results. 
Averaged  mass  concentration  profiles  across  the  center  of  the  vortex  as 
well  as  through  the  laboratory  fireball  wake  are  discussed.  An  iso-concen¬ 
tration  picture  which  provides  the  complete  averaged  spatial  distribution 
of  the  low  density  gas  in  a  turbulent  laboratory  event  is  presented.  Much 
of  the  data  acquired  during  the  course  of  this  study  is  discussed,  but 
for  completeness,  all  acquired  data  was  tabulated  and  is  included  in 
Appendix  E. 
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2.  EXPERIMENTAL  FACILITY 


A  test  facility  consisting  of  several  primary  components  including 
a  high  pressure  test  tank  with  optical  ports,  equipment  for  the  remote 
formation,  sensing,  and  controlled  bursting  of  soap  bubbles  filled  with  low 
density  gas,  and  the  necessary  diagnostic  equipment  to  obtain  density 
field  measurements  and  flow  visualization  photographs  of  the  internal  and 
external  vortex  flow  were  designed,  developed,  and  utilized  in  the 
performance  of  this  experimental  study. 

2.1  HIGH  PRESSURE  TANK  AND  AUXILIARY  EQUIPMENT 

The  high  pressure  tank  is  an  ASME  pressure  vessel  with  a  diameter 
of  54  inches  and  a  length  of  14  feet.  A  photograph  of  the  test  tank 
installed  in  the  TRW  Fluid  Mechanics  Laboratory  is  shown  in  Figure  1.  This 
facility  was  designed  to  operate  in  a  pressure  range  from  vacuum  conditions 
to  10  atmospheres. 

A  system  was  designed  and  developed  to  renotely  form  soap  bubbles 
filled  with  low  density  gas  (relative  to  the  ambient  higher  density  gas) 
in  the  high  pressure  test  tank.  Helium  was  used  as  this  low  density  gas 
when  air  was  used  as  an  ambient  medium,  while  nitrogen-filled  bubbles  were 
used  in  the  higher  Reynolds  number  SF6  tests.  Positive  displacement  pistons 
and  a  valving  system  were  used  to  meter  the  amount  of  low  density  gas  and 
soap  solution  to  the  bubble  tube  to  accurately  control  the  bubble  size. 

A  small  electrode  spaced  about  20  mils  from  the  bubble  tube  with  a 
200  microfarad  capacitor  discharge  system  was  initially  used  to  burst  the 
bubbles.  As  a  result  of  early  asymmetries  present  in  vortices  generated 
by  this  initial  electrical  discharge  bubble  breaking  system,  a  mechanical 
breaking  technique  was  developed  which  repeatedly  produced  a  symmetric 
helium  release.  This  technique  consisted  of  remotely  dropping  a  125  mil 
stainless  steel  sphere  from  an  electromagnet  mounted  inside  the  tank 
through  the  top  center  of  the  bubble.  Tests  with  no  bubble  have  shown  that 
no  more  -han  1/16"  dispersion  from  the  target  ever  existed  and  interferograms 
taken  just  milliseconds  after  arrival  of  the  sphere  at  all  test  pressures 
indicated  that  this  system  provided  a  symmetric  helium  release. 


9 


Figure  1.  High  Pressure  Test  Facility,  Holocamera,  and  Support  Equipment 


The  bubble  formation  head  and  electrode  breaker/sensor  were  designed 
as  an  integral  unit  so  they  could  be  traversed  vertically  within  the  tank 
through  a  1-inch  diameter  access  flange  and  seal.  The  bubble  tube  position 
could  be  varied  from  the  center  to  30  inches  below  the  horizontal  viewing 
ports.  As  a  result  of  the  complexity  of  the  optical  instrumentation,  and 
in  particular  the  interferometry  apparatus,  it  was  essential  that  they 
remained  fixed  relative  to  the  tank  viewing  ports  which  necessitated  the 
vertical  relocation  of  the  event  initialization  point  in  order  to  provide 
data  at  the  various  desired  rise  positions. 

2.2  DIAGNOSTIC  EQUIPMENT 

In  order  to  completely  specify  the  fluid  state  and  motion,  the 
following  must  be  known  for  a  non-reacting  gas:  two  thermodynamic  variables, 
such  as  pressure  and  temperature;  a  variable  of  the  motion,  such  as  velocity, 
and  the  composition  of  1  he  field.  Measurements  of  the  tank  pressure  and 
temperature,  outside  the  induced  flow,  specified  the  pressure  and  temperature 
of  the  field  since  the  induced  fluid  motion  was  incompressible  (i.e.,  low 
velocity  and  adiabatic).  Of  the  remaining  two  unknowns,  the  density  or 
composition  field  was  measured  optically  using  holographic  interferometry 
and  although  the  velocity  field  was  not  measured  directly,  an  optical 
particle  tracking  technique  was  used  to  provide  information  on  the  internal 
unsteady  flow  structure.  In  addition,  a  shadowgraph  system  was  used  to 
provide  data  on  the  external  shape  and  position  of  the  events  as  a  function 
of  time.  Optical  techniques  were  adopted  for  this  study  because  of  three 
primary  difficulties  involved  with  flowfield  probes.  First,  interpretation 
of  data  obtained  from  probes  in  a  highly  three-dimensional  flow  is  almost 
impossible.  Also,  probe  interference  effects  are  sizable  in  an  incompressible 
flow  such  as  this,  and  finally  the  cost  and  time  consumed  in  reconstructing 
a  time  dependent  flowfield  from  local  measurements  was  prohibitive. 
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2-2.1  Holographic  Interferometer 

Pulsed  laser  holographic  interferometry  was  used  to  measure  the 
mean  density  field  of  the  rising  toroids.  The  laser  holograph*  used  in 
this  program  was  similar  in  design  to  one  used  in  previous  fluid  mechanic 
studies  at  TRW  but  was  considerably  larger  in  size  to  accommodate  the 
large  high  pressure  tank.  The  main  feature  of  the  holocamera  is  that  it 
records  orthogonal  views  simultaneously  across  the  tank  test  section. 
Consequently,  this  system  is  capable  of  recording  any  transmission  subject 
in  an  8-inch  diameter  area  in  line  with. the  4-  by  5-inch  photographic 
plates  and  the  large  18-inch  objective  lens  at  the  opposite  ports.  The 
advantages  of  using  this  holographic  interferometer  as  compared  to  the 
standard  Mach  Zehnder  interferometer  include  the  three-dimensional  viewing 
and  depth  focusing  properties  of  holographic  interferograms,  and  thi  fact 
that  precision  optics  and  accurate  alignment  are  not  required  with 
holographic  interferometry  since  a  common  optical  path  is  used  fov'  both  the 
scene  and  comparison  beams.  A  schematic  of  the  complete  holocamera 
identifying  all  primary  optical  components  is  presented  in  Figure  2.  A 
photograph  of  the  holocamera  together  with  the  test  tank  was  shown  in 
Figure  1 . 

In  this  holocamera  the  primary  beam  from  the  laser  is  split  into  a 
reference  and  a  scene  beam.  The  scene  beam  is  directed  through  the  event 
to  the  plates  whereas  the  reference  beam  is  directed  on  the  plates  without 
passing  through  the  disturbance.  The  prism  plate  and  lenses  shown  in  Figure  2 
are  used  to  spatially  and  temporally  match  the  scene  and  reference  beams 
over  the  plates.  A  Q-switc:hed  solid  state  ruby  laser  with  an  energy  of 
about  1  joule  and  pulse  length  of  about  100  nanoseconds  was  used  to  expose 
both  the  vertical  and  horizontal  beam  holograms  which  were  recorded  on 
Agfa-Gavaert  8E75  plates. 


*The  term  laser  holograph  pertains  to  both  the  laser  cavity  or  oscillator, 
and  all  other  optical  components  which  comprise  the  holocamera. 
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Figure  2.  Schematic  Diagram  of  Holocamera 


Interferograms  were  obtained  in  this  common  path  interferometer  by 

sequentially  recording  the  comparison  and  scene  beams  on  the  same  plates. 

The  double  exposure  of  the  plates  by  the  two  pulses,  the  one  of  the 

undisturbed  volume  and  the  other  of  the  same  volume  with  the  test  scene 

introduced,  produced  the  interferograms.  The  principle  of  operation  of 
any  interferometer  is  that  optical  interference  occurs  between  two  coherent 
beams  passing  through  media  of  different  density.  When  optical  paths 
differ  by  an  integral  multiple  of  the  wavelength  of  the  monochromatic  light 
source,  X,  (AL  =  nX,  n  =  0,  1,  2.  3  ...)  the  beams  reinforce  one  another. 
Conversely,  when  the  optical  paths  differ  by  odd  multiples  of  X/2  (aL  =  x/2, 
3X/2,  5x/  2  ...),  complete  interference  occurs.  Reinforcement  shows  upas 
light  fringes  on  the  holographic  plate  whereas  complete  interference  creates 
dark  fringes. 

When  the  finite  fringe  technique  is  used,  the  plate  holder  and  its 
adjacent  mirror  are  rotated  between  event  and  comparison  exposures  in  such 
a  manner  that  the  optical  paths  cf  the  scene  and  reference  beams  change 
uniformly  over  the  entire  plane  of  the  holographic  plate,  creating  a  raster 
of  parallel  fringes.  The  spacing  of  these  fringes  is  dependent  on  the 
amount  of  rotation  of  the  plate  holder  and  mirror.  The  more  the  mirror  and 
plate  holder  are  displaced,  the  greater  is  the  path  mismatch,  and  consequently, 
the  interference  fringes  become  narrower.  Whatever  additional  shift  is 
caused  by  the  test  disturbance  shows  up  as  fringes  displaced  from  the  uniform 
raster.  When  the  infinite  fringe  technique  is  implemented,  no  optical 
components  are  rotated  between  exposures.  Only  the  disturbance 
in  a  blank  background  appears,  hence  the  name  infinite  fringe.  Various 
examples  of  finite  and  infinite  fringe  interferograms  were  presented  in 
References  1  and  2  and  appear  in  this  report  in  Figures  3  through  7.  A 
further  discussion  of  the  acquisition  and  interpretation  of  both  finite 
and  infinite  fringe  interferograms  and  the  advantages  and  disadvantages 
of  each  is  included  in  Section  3. 
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'  2.2.2  Particle  Tracking  Apparatus 

Equipment  providing  the  capability  of  obtaining  multiple  exposure 
photographs  of  a  rising  vortex  with  its  initial  bubble  seeded  with  a 
particulate  tracer  was  developed  and  used  to  provide  data  on  the  internal 
vortex  structure. 

2. 2. 2.1  Particulate-Seeded  Bubbles 

Several  different  types  of  flow  tracers  were  tested,  such  as  talcum 
powder,  fine  aluminum  flakes,  smoke,  and  hollow  epoxy  microspheres.  Smoke 
was  discarded  since  it  could  not  be  generated  at  the  high  operating  pressures, 
talcum  powder  could  not  be  photographed,  and  aluminum  flakes  prohibited 
repeatable  bubble  formation.  Consequently,  epoxy  microspheres  from  10  to 
250  microns  in  diameter  were  incorporated  as  the  tracing  particulate.  It 
was  found  that  these  microspheres  were  both  a  good  tracer  and  were  sufficiently 
visible  when  illuminated  to  be  photographed.  An  analysis  presented  in 
Appendix  A  confirms  the  effectiveness  of  these  particles  as  fluid  tracers. 

2. 2. 2. 2  Camera/Flash  System 

A  35  mm  Nikon  F  with  motor  drive  having  an  f5.6  macro  lens  was 
selected  as  the  optimum  available  camera  to  acquire  the  particle  tracking 
photographs .  Kodak  Tri-X  Pan  (ASA  400)  was  chosen  for  film,  and  a  special 
processing  technique  using  Acufine  developer  was  incorporated  which  effec¬ 
tively  increased  the  speed  of  the  film  to  ASA  1250,  noticeably  enhancing 
the  resolution  above  that  obtainable  using  standard  developers. 

A  7.5kV  power  supply  operating  at  an  energy  level  of  about  50  joules 
was  developed  to  rapidly  charge  a  2.5yF  capacitor  bank  which  when  discharged, 
ignited  an  EG&G  FX98C-3  (3"  length,  70mm  diameter)  flash  lamp.  The 
capacitor  discharge  was  triggered  by  a  standard  square  wave  signal  generator 
at  the  desired  frequency.  An  oscilloscope  sweep  triggered  by  the  release 
of  the  bubble  breaking  sphere  from  the  electromagnet  facilitated  accurate 
monitoring  of  the  flash,  enabling  the  time  at  which  each  exposure  was  taken 
after  bubble  breaking  to  be  accurately  known.  The  flash  was  masked  except 
for  two  aligned  slits  to  provide  a  pulsed  vertical  sheet  of  light  directly 
above  the  bubble  tube  and  normal  to  the  line  of  sight  of  the  camera.  This 
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illumination  system  permitted  particulate  to  be  visible  in  only  a  cross- 
sectional  slice  (less  than  1/2"  wide)  through  each  event  allowing  detailed 
visualization  of  the  internal  flow  structure. 

2.2.3  Shadowgraph  System 

In  general,  any  light  ray  passing  through  a  region  whose  index  of 
refraction  varies  is  displaced,  deflected,  and  changed  in  phase  relative  to 
the  path  it  would  have  taken  if  the  index  of  refraction  were  to  have  been 
constant.  If  a  collimated  light  boam  is  passed  through  a  disturbed  region, 
the  image  one  sees  on  a  screen  on  the  other  side  of  the  light  source  from 
the  disturbance  becomes  a  shadowgraph.  Using  the  approximation  of  small 
deflections,  for  most  gases  it  may  be  shown  that  the  change  in  light 
intensity  on  the  screen  due  to  the  disturbance  is  proportional  to  the 
second  derivative  of  the  density  within  the  disturbed  region  (see  Reference  3 
for  example). 

The  shadowgraph  system  used  in  this  study  utilized  a  point  light 
source  (500  watt  Mercury  arc  lamp)  to  the  right  of  the  large  lens  adjacent 
to  the  tank  window  in  Figure  2.  The  lens  collimated  the  beam  and  a  trans¬ 
lucent  vellum  screen  was  placed  on  the  left  tank  window.  A  16  mm  Bell  and 
Howell  movie  camera  with  a  wide-angle  lens  using  Eastman  4-X  movie  film 
(ASA  500)  at  24  frames/second  was  used  to  record  the  events. 
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3.  DATA  ACQUISITION  AND  INTERPRETATION 


The  purpose  of  this  section  is  to  discuss  both  the  procedure  for 
obtaining  infinite  fringe  and  finite  fringe  holographic  interferograms  and 
the  physical  interpretation  of  these  forms  of  data.  As  well,  the  acquisition 
and  interpretation  of  particle  tracking  photographs  and  the  acquisition  of 
shadowgraph  data  are  discussed. 


3.1  ACQUISITION  OF  INTERFEROMETRY  DATA 

The  procedure  which  was  used  to  obtain  a  set  of  infinite  fringe 
interferograms  of  a  simulated  fireball  consisted  of  the  following  primary 
steps : 

1)  Plate  holders  containing  the  special  fine  grain,  low  speed 
photographic  plates  were  loaded  in  the  horizontal  and 
vertical  plate  holder  receptacles. 

2)  The  slides  covering  the  plates  were  drawn  after  the  laser  and 
bubble  blowing  apparatus  had  been  checked  and  all  sources  of 
illumination  in  the  test  area  had  been  eliminated. 

3)  The  laser  power  supply  was  charged  and  a  soap  bubble  filled 
with  low  density  gas  was  blown  to  a  predetermined  size. 

4)  About  five  seconds  after  its  formation,  the  bubble  was  burst. 

5)  Simultaneous  with  the  bubble  bursting,  an  electrical  signal, 
first  sent  through  a  time  delay  circuit,  activated  the  laser 
firing  mechanism  which  recorded  the  test  scene  on  the  plates. 
The  time  delay  was  adjusted  to  allow  the  helium  cloud  to  rise 
a  predetermined  amount. 

6)  A  comparison  scene  was  recorded  on  the  plates  by  firing  the 
laser  a  second  time  about  30  seconds  after  the  event  sce.ie. 
This  delay  ensured  that  the  disturbance  was  well  outside  of 
the  viewing  region  resulting  in  a  uniform,  constant  density 
comparison  scene. 

7)  The  photographic  plates  were  removed  and  developed  using 
special  photographic  developing  agents. 

8)  The  holograms  were  reconstructed,  i.e.,  observed  and 
photographed  by  illuminating  the  hologrars  with  a  CW  laser. 

9)  Large  scale  working  photographs  were  then  produced  for  use 
in  the  data  reduction  phase  as  described  in  Section  4.1. 
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A  horizontal  view  finite  fringe  interferogram  was  obtained  by  following 
the  exact  procedure  as  outlined  above  except  for  the  addition  of  one  step. 

This  consisted  of  rotating  the  horizontal  beam  plate  holder  and  adjacent 
mirror  a  small  precise  amount  between  steps  5  and  6  above,  i.e.,  between 
recording  the  test  and  the  comparison  scenes  on  the  plate.  The  amount  of 
this  rotation  determines  the  spacing  of  the  undisturbed  fringes. 

3.2  INTERPRETATION  OF  HOLOGRAPHIC  INTERFEROGRAMS 

Interpretation  of  both  infinite  fringe  and  finite  fringe  interfero- 
grams  are  included  in  this  section  since  both  types  o,:  data  were  acquired 
and  analyzed  during  the  course  of  this  study. 

3.2.1  Infinite  Fringe  Interferograms 

An  example  of  an  infinite  fringe  interferogram  is  presented  in 
Figure  3.  This  horizontal  view  hologram  was  obtained  with  an  unbroken 
one-inch  diameter  helium  bubble  at  a  tank  pressure  of  one  atmosphere.  Each 
fringe  in  this  type  of  interferogram  corresponds  to  a  locus  of  constant 
optical  path  so  interpretation  for  data  reduction  purposes  consists  of 
numbering  and  locating  the  fringes  along  any  plane  of  interest  through  the 
fireball.  These  fringe  numbers  and  locations  can  directly  be  used  to 
determine  the  integrated  density  or  helium  concentration  change  through  the 
fireball  along  that  optical  path.  Details  of  the  computational  procedure 
required  to  obtain  concentration  profiles  from  these  fringe  number  data  are 
presented  in  Section  4.1.  In  the  sample  infinite  fringe  interferogram 
presented,  the  density  is  constant  inside  the  bubble  and  lower  than  the 
constant  density  region  outside  the  bubble  by  about  a  factor  of  7.  Consequently, 
the  interference  fringes  which  consist  of  concentric  circles  with  decreasing 
spacing  as  their  radii  increase  depict  only  changes  in  optical  path  length 
through  the  helium  sphere. 

Another  example  of  an  infin'te  fringe  interferogram  is  shown  in 
Figure  4.  This  interferogram  shows  a  laminar  vortex  formed  from  an  initial 
one-inch  diameter  helium  bubble  at  one  atmosphere  pressure  after  a  rise  of 
about  4  initial  bubble  diameters.  The  fringe  pattern  in  this  figure  is  more 
typical  of  that  expected  in  this  study  in  that  the  fringe  geometry  and  spacing 
result  from  a  toroidal  cloud  containing  a  mixture  of  helium  and  air. 
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nite  Fringe  Interogram 


Figure  5.  Typical  Turbulent  Infinite  Fringe  Interferog 


A  turbulent  vortex,  typical  of  the  higher  Reynolds  number  results  in 
this  study,  is  shown  in  Figure  5  as  a  third  example  of  an  infinite  fringe 
inter ferogram.  This  event  was  formed  from  an  initial  one-inch  diameter  helium 
bubble  at  eight  atmospheres  pressure  after  a  rise  of  about  5  initial  bubble 
diameters  (Reynolds  number  a  factor  of  8  larger  than  the  event  in  Figure  4). 

It  is  obvious  by  comparing  this  event  with  the  lower  Reynolds  number  hologram 
in  Figure  4  that  the  higher  Reynolds  number  fireball  contains  small  scale 
turbulent  structure,  not  esent  in  the  laminar  cast. 

3.2.2  Finite  Fringe  Interferograms 

An  example  of  a  finite  fringe  interferogram  is  presented  in  Figure  6. 
This  interferogram  was  obtained  at  the  same  test  conditions  as  the  infinite 
fringe  interferogram  in  Figure  4  but  at  a  higher  rise  position  (6  initial 
bubble  diameters).  In  this  type  of  hologram,  the  test  disturbance  causes  a 
displacement  or  shift  in  the  fringes  relative  to  the  undisturbed  portion 
of  the  test  section.  An  increase  in  optical  path  is  shifted  in  the  opposite 
direction  to  a  decrease  in  optical  path.  Consequently,  no  ambiguity  can 
occur  in  regard  to  the  fringe  shift  direction.  The  magnitude  of  the  fringe 
shift  in  this  type  of  hologram  can  be  directly  related  to  fringe  number  in 
an  infinite  fringe  hologram  and  therefore  to  an  integrated  change  in  density 
or  concentration  along  the  optical  path. 

A  second  example  of  a  finite  fringe  interferogram  is  presented  in 
Figure  7.  This  hologram  of  a  turbulent  vortex  was  obtained  at  the  same 
conditions  as  the  infinite  fringe  interferogram  of  Figure  5  but  at  a  slightly 
higher  rise  position  of  6  initial  bubble  diameters.  As  demonstrated  by  this 
example,  finite  fringe  holograms  with  easier-to-reduce  continuous  fringes 
can  be  obtained  in  small  events  with  high  helium  concentrations  by  increasing 
the  amount  of  mirror  and  plate  rotation  between  recording  the  test  scene 
and  comparison  scene  on  the  hologram  plates. 

3.2.3  Advantages/Disadvantages  of  Infinite  Fringe/Finite  Fringe 

The  advantages  of  the  infinite  fringe  arrangement  are  twofold.  First, 
the  overall  setup  and  care  of  the  optics  is  much  simpler  since  path  matching 
of  the  two  beams  is  not  altogether  critical.  Secondly,  the  infinite 


Figure  6.  Typical  Laminar  Finite  Fringe  Interferogram 


Figure  7.  Typical  Turbulent  Finite  Fring 


Interferogram 


fringe  interferogram  provides  an  excellent  outline  of  the  test  event  and  a 
clear  physical  picture  of  the  internal  flowfield. 

However,  the  advantages  of  the  finite  fringe  technique,  especially  in 
providing  a  wide  range  of  test  conditions,  outweigh  for  the  most  part  the 
merits  of  the  infinite  fringe  technique.  When  there  exist  relatively  large 
density  gradients  in  the  flow  (i.e.,  very  many  fringes  of  shift  along  a 
small  length),  and  when  the  flow  is  turbulent,  the  fringes  tend  to  be 
very  irregular  in  shape  and  difficult  to  delireatt  in  an  infinite  fringe 
iocerferogram.  Also,  it  is  not  known  internally  whether  a  fringe  shift 
means  an  increase  or  decrease  in  density.  Consequently,  an  a  priori 
knowledge  of  the  flowfield  is  a  prerequisite  to  use  that  technioue.  Even 
though  the  infinite  fringe  approach  is  more  simple  and  provides  better 
visual  insight  to  the  flow,  the  flexibility  of  the  finite  fringe  technique 
(rotating  elements  in  an  amount  commensurate  with  the  character  of  the  flow) 
proved  \.c  be  most  advantageous.  Without  its  implementation,  the  wide  range 
of  Reynolds  number  variation  at  identical  stages  of  vortex  development  and 
especially  the  acquisition  of  high  Reynolds  number  data  could  not  have  been 
accomplished. 
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4.  DATA  REDUCTION 


The  objective  of  this  section  is  to  describe  in  detail  techniques 
which  were  used  in  reducing  nolographic  interferometry  data  (both  infinite 
fringe  and  finite  fringe)  and  shadowgraph  data  as  they  applied  to  this  study. 

4.1  HOLOGRAPHIC  INTERFEROMETRY  DATA 

The  process  of  achieving  the  fina.l  result,  namely  determining  the 
mean  density  and  mass  concentration  fields  from  the  holographic  records 
consisted  of  two  primary  steps.  The  first  step  consisted  of  fringe 
interpretation,  determination,  and  averaging.  Once  an  average  fringe 
distribution  was  determined  for  a  large  set  of  repeat  data,  computer  reduction 
of  the  density  and  mass  concentration  fields  was  able  to  be  accomplished.  This  basic 
approach  has  been  used  in  a  number  of  former  studies  at  TRW  (References  4, 

5  and  6)  but  the  averaging  technique  is  unique  to  this  study. 

The  fringe  shift  S  measured  from  the  interferograms  is  related  to 
the  index  of  refraction  n  (from  which  density  and  mass  concentration  are 
determined)  according  to  the  following  integral  equation: 

s  =  !/  (n  -  n»)ds 
-oo 


where  X  =  wavelength  of  light  in  vacuum 

ds  =  differential  path  length  along  light  ray 
and  subscript  °°  corresponds  to  undisturbed  conditions. 

The  problem  reduces  to  solving  this  integral  equation  for  the  change 
in  index  oi  refraction,  (n  -  nj,  as  a  function  of  fringe  shift,  S.  The 
mechanics  of  inverting  the  integral  will  be  discussed  in  detail  in  Appendix 
B. 

This  solution  is  accomplished  by  restricting  consideration  to  an 
ax i symmetric  geometry.  Naturally,  to  use  this  technique  to  produce  results 
that  are  meaningful,  the  input  fringe  data  must  originate  from  an  axisymmetric 
event.  To  obtain  an  axisymmetric  event,  some  type  of  average  of  repeat 
sets  of  fringe  distributions  must  be  obtained.  One  could  take  an  average 
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in  such  a  way  to  yield  a  probabil  ity-type  result  for  density  at  a  specified 
location  in  space,  i.e.,  a  tank-fixfd  approach.  This  would  result  in  a  very 
smeared  out*  vortex  which  in  no  way  would  resemble  a  typical  vortex. 

Instead,  an  approach  was  developed  whereby  a  coordinate  system  was  established 
for  each  individual  vortex  based  on  the  location  of  three  points  common  to 
all  vortices,  but  located  in  general  differently  from  event  to  event.  This 
approach  retained  the  important  physical  features  of  the  flow  and  produced 
a  symmetric  vortex  which  was  a  representative  average  vortex  rather  than  a 
smeared  out  vortex  average.  This  averaged  fringe  distribution  forms  the 
topic  for  the  next  section. 

4.1.1  Averaging  Radial  Fringe  Shift  Distributions 

To  outain  an  averaged  radial  fringe  distribution  for  a  particular  set 
of  repeat  interferograms,  several  basic  steps  were  undertaken.  First,  coor¬ 
dinates  were  established  on  the  interferograms,  and  individual  fringe 
profiles  were  constructed.  The  data  were  then  all  collapsed  about  points 
characteristic  to  all  vortices  in  general.  Once  this  had  been  accomplished, 
an  average  S(r)  distribution  was  calculated,  and  being  satisfactorily 
symmetric,  was  then  used  as  the  input  to  the  density  field  and  mass  concentra¬ 
tion  field  calculations. 

4.1 .1 .1  Establishment  of  Coordinates  on  Individual  Interferograms 

An  infinite  fringe  interferogram  was  visually  examined  to  determine 
the  position  of  the  saddle  point  by  counting  fringes  inward  from  the  edge  of 
the  vortex.  In  this  connotation,  the  saddle  point  is  defined  as  the  point 
where  the  fringe  number  (see  Section  4. 1.1. 2)  simultaneously  attains  a 
maximum  value  in  an  axial  direction  and  a  minimum  value  in  the  radial 
direction.  Two  fringe  maxima  positions  (generally  having  different  fringe 
numbers)  were  determined  in  the  same  manner  and  an  x-axis  was  drawn  through 
them.  The  y-axis  was  then  drawn  orthogonal  to  the  x-axis  and  through  the 
saddle  point.  The  point  where  the  x-  and  y-axes  intersected  defined  the 

*In  the  sense  that  the  averaged  event  would  be  much  larger  than  any  individual 
vortex  and  the  important  edge  gradients  would  be  reduced. 
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vortex  coordinate  origin,  or  for  the  sake  of  brevity,  simply  the  origin. 
Reference  to  Figure  8  shows  that: 


xL£  =  distance  from  origin  to  first  dark  fringe  on  left  side 

XLM  =  distance  from  origin  to  fringe  maximum  on  left  side 

x^  =  distance  from  origin  to  fringe  maximum  on  right  side 


x 


RE 


=  distance  from  origin  to  first  dark  fringe  on  right  side 


distance  from  origin  to  saddle  point 


Tank-fixed  parameters  of  purely  statistical  interest  are: 


n  =  rotation  of  x-axis  from  true  horizontal 


zQ  =  vertical  rise  distance  of  origin  from  top  of  bubble  tube 

x„  =  horizontal  distance  from  origin  to  center  of  bubble  tube 
o  3 


Figure  9  represents  a  schematic  diagram  of  a  finite  fringe  inter- 
ferogram  which  illustrates  the  placement  of  coordinates.  Note  that  locating 
the  maxima  and  saddle  point  was  not  accomplished  by  counting  fringes  inward 
from  the  edges,  but  by  determining  relative  displacements  of  fringes  from 
their  undisturbed  positions.  These  features  do  not  become  immediately 
obvious  as  was  the  case  with  finite  fringe  interferograms,  and  to  identify 
them  correctly  requires  some  hunting.  Once  these  three  points  have  been 
correctly  determined,  all  coordinates  become  established  in  an  identical 
fashion  to  the  infinite  fringe  coordinates. 

One  set  of  finite  fringe  interferograms  in  particular  was  more 
thoroughly  analyzed  for  the  purpose  of  constructing  a  map  of  the  entire  vortex 
vortex  field.  For  this  set  (10  atmospheres  pressure,  N2  bubbles 
in  an  SF6  environment,  10.4  initial  bubble  diameters  rise,  and  an  initial 
bubble  diamet°r  of  1/2  inch),  four  additional  cuts  through  the  vortex  as 
well  as  three  cuts  through  the  fireball  wake  were  made.  The  vertical 
positions  of  the  four  additional  vortex  cuts  parallel  to  the  central  cut 
were  normalized  with  respect  to  where  the  top  of  the  vortex  intersected  the 
y-axis,  y^,  and  where  the  bottom  of  the  vortex  intersected  the  y-axis,  yB> 
These  cuts  were  made  at  1/3  y^,  2/3  yt,  1/3  yg,  and  2/3  yB>  For  the  1/3  yt 
and  1/3  yB  cuts,  the  two  maxima  and  central  minimum  were  identified  in  an 
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Figure  8.  Schematic  of  Interferometry  Data  Reduction 
Infinite  Fringe 


identical  fashion  to  the  central  cut.  For  the  2/3  yt  and  2/3  yg  cuts,  more 
often  than  not  maxima  could  not  be  detected  and  consequently  the  central 
minimum  became  the  only  criterion  for  coordinate  selection.  No  bias  from  the 
central  cut  was  placed  on  locating  the  central  minimum  or  the  maxima  to 
again  preserve  the  physical  flow  features.  The  motivation  behind  allowing 
these  coordinates  to  float  was  that  in  the  mean,  averaged  over  a  sufficient 
number  of  events,  they  would  coincide  with  those  of  the  central  cut  (results 
show  that  they  did  indeed). 

4 . 1 . 1 . 2  Obtaining  Fringe  Profiles  from  Individual  Interferograms 

For  infinite  fringe  interferograms,  the  center  of  each  fringe  was 
located  and  numbered  where  it  intersected  the  x-axis  defined  in  the  previous 
section.  The  numbering  started  with  S  =  1/2  for  the  outermost  dark  fringe 
(corresponding  to  the  half-wavelength  shift  in  optical  path  discussed 
earlier),  S  =  1  for  the  next  inward  light  fringe,  and  so  on,  increasing  by 
1/2  each  time  a  new  fringe  was  encountered.  Care  had  to  be  taken  to  give  the 
same  fringe  the  same  number  when  it  intersected  the  x-axis  more  than  once. 

What  resulted  were  m  ordered  pairs  (x^,  Sij.)  j  =  1 ,  2,  . . .  m  for  every  ith 
interferogram,  i  =  1,  2,  ...  n. 

Because  of  the  nature  of  the  finite  fringe  interferogram,  this  step 
in  the  data  reduction  procedure  was  considerably  different  than  the  above 
for  infinite  fringe  data.  Reference  again  to  Figure  9  shows  that  all 
fringes  were  first  identified  by  number  (not  to  be  confused  with  S  =  1/2,  1, 

1  1/2,  etc.  for  infinite)  in  the  undisturbed  regions  both  to  the  left  and 
to  the  right  of  the  event.  A  mark  was  made  wherever  the  center  of  any 
fringe  intersected  the  x-axis,  and  each  mark  was  identified  by  its  fringe 
identification  number  by  following  its  fringe  into  the  undisturbed  region 
where  it  was  numbered.  Then,  the  distance  normal  to  the  parallel  fringes 
from  the  disturbed  fringe  to  where  its  undisturbed  fringe  would  be  if  there 
were  no  disturbance,  d..,  was  measured,  and  divided  by  the  fringe  spacing, 

*  J 

AL.,  to  give  S  in  terms  of  number  of  fringes  of  shift.  This  S  was  completely 
analogous  to  the  S  obtained  by  counting  fringes  inward  from  the  edge  for 
infinite  fringe  interferograms. 

It  should  be  noted  that  all  that  follows  in  the  remainder  of  Section  4.1 
applies  equally  to  both  infinite  ana  finite  fringe  data. 
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4. 1.1. 3  Data  Averaging 

In  the  process  of  establishing  an  average  vortex,  i  i  was  essential 
to  preserve  the  important  physical  features  by  normalizing  the  characteristic 
dimensions  of  each  individual  vortex.  The  dimensions  selected  for  this 
purpose  were  the  locations  of  the  fringe  maxima  and  the  left  and  right  edges 
relative  to  the  origin. 

To  transform  each  event  into  a  system  of  coordinates  whereby  maxima 
collapsed  onto  the  mean  maxima  and  edges  collapsed  onto  the  mean  edges, 
the  general  linear  transformation  x. ,  =  ax..  +  b  was  used.  First,  the 

1  J  1  J 

average  locations  of  the  edges  and  the  maxima  were  computed  using 

n 

*LE  '  «  £  XLEi 
1  n 

XLM  =  n  S  XLM. 
i=l  1 

n 

XRM  =  n  S  XRM. 
i=l  1 

^  n 

XRE  =  n  S  XRE . 

i=l  i 

The  event  was  then  divided  into  four  regions  as  shown  below  and 
the  normalization  was  accomplished  separately  within  each  region. 


I 

I 

1  1 

1 

1 

1  1 

l  Reqion 

1  !  Region  2 

Region  3  [  Region  4  ’ 

XLE . 

i 

xLMj 

1  ( 

\  XRM-f  xREi 

The  constants  for  a  and  b  in  the  general  linear  transformation  were  determined 
by  matching  the  endpoints  within  each  region  to  the  averaged  endpoints. 

For  example,  in  region  1,  at 


x 


ij 


and,  at  x 


ij 
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iL  4-  U 

For  the  jtn  data  point  on  the  iLn  interferogram,  the  new  coordinates 


x..  were: 

I  J 


In  region  1, 


XLE.  -  xij  -  xLMi 


x  • .  = 


XLE  XLM  +  -  XLE  "  XLM 

A  •  •  '  A|  u  ”  ..  A| 


'ij  X,  r  -  x,  M  ‘ij  LM  X,  r  -  X,  M  LM. 


!Ei  ALM. 


kLE.  LM.  i 


In  region  2, 


*LH.  i  x,j  i.  0 


7..  -  H-JL*.. 

^  XLM.  ^ 


In  region  3, 


0  <  xi-j  <  xrM. 


X..  =  -M_x 
1J  XRM .  1J 


In  region  4, 


XRMi  -  xij  -  XRE. 


x..  -  xr£  :  xr»  +  x 

1J'  XRE(  "  "rM^  '•>  RM 


XRE  "  XRM 
XREi  "  xRMi  XRMi 


In  cuts  where  normalization  was  accomplished  only  with  respect  to 
the  edges  e.g.,  (2/3  y^,  2/3  yB,  wake  cuts,  etc.),  the  transformations  were 
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..  .  - 


for  the  left  side, 


ij 


LE 

<LEi 


Xij 


and  for  the  right  side, 


ij 


RE 

<REi 


Xij 


Once  all  (x..,  S..)  had  been  transformed  to  (x..,  S..),  there  existed  n 
_ 1  j  I  j  i  j  __  ;  j 

S1d  vs.  x^.  profiles  with  all  edges  occurring  at  xLE  and  xRE,  all 
maxima  at  and  XR^,  and  the  central  minimum  at  x  =  0.  These  profiles  had 
to  be  averaged  with  respect  to  S  at  constant  x  to  yield  a  mean  S  vs.  x 
profile.  Since  every  i  n  distribution  did  not  have  a  discrete  value  of  S  at 
the  same  x,  convenient  xk  were  chosen  to  be  used  for  each  ith  event  (say 

xk  =  U  XLE’  ^  =  2»  ...  L)  and  linear  interpolation  for  Sik  was  used 

between  (x, j .  Sfj)  and  (xjj+].  S.j+1). 


S1k  '  (7k  -  V  +  Sij 

ij+1  "  xi j 


The  input  fringe  numbers  to  the  Abel  inversion  data^reduction  routine  for 

an  axisymmetric  average  vortex  were  simply  S.  =  -  J2  S ..  for  k  =  1 ,  2,  .. 

k  n  j  _  j  ik 

The  S  vs.  x  for  both  the  left  and  right  sides  were  then  computed.  For  all 
practical  purposes,  the  two  sides  coincided  (see  Figure  17  for  a  typical 
example)  when  the  number  of  events  averaged  was  sufficiently  large. 


L. 


4.1.2  Calculation  of  Density  and  Mass  Concentration  Fields 

Given  the  fringe  distribution  for  an  axisymmetric  event,  it  is  a 
straightforward  matter  to  calculate  a  density  and  mass  concentration  dis¬ 
tribution  along  the  same  cut. 
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4. 1.2.1  Abel  Inversion  of  Fringe  Shift  Equation 

Recall  the  general  equation  for  fringe  shift 


S  =  1  /  -  "Jds 

•'-CO 

Assuming  the  axi symmetric  geometry  below 


it  is  seen  that  the  incremental  path  length  ds  becomes  d[/r2  -  y2],  and 
since  n  =  n^  for  all  r  >  rQ  and  that  the  event  consists  of  two  equal 
halves,  the  fringe  shift  equation  for  an  axi symmetric  event  becomes 


(1) 


To  solve  for  n  -  n^,  one  notes  that  this  equation  is  the  same  as  Vol terra's 
integral  equation  of  the  first  kind,  with  the  exception  that  the  limits  of 
integration  are  different.  The  details  of  the  solution  (Reference  7)  are 
rather  length  and  are  presented  in  Appendix  B.  It  becomes 


(n  -  "Jr.y  =  -£/r° 

w  ^  w 


dS  j 
3r  dr 


7y  (r2  -  y2)1/2 
Break  this  up  into  n  integrals  over  n  intervals 


(2) 


n  -  n  ) 


r=y 


iii  /• 

1  k=l  r. 


n-l  -rk+l  <3F>  dr 


(r2  -y?) 


1/2 


i  =  0,  1,  2,  3 


Assume  that  in  each  kth  zone  (^)  is  constant.  Then, 
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n  -  nJr=y  =  /  k  T7 — r-2"  i/2 

k=i  k  (r2  -  r?)  ' 


n-1 


-k_ 

TT 


k+1 


-  S, 


R  rk+l  ‘  rk 


log 


r  +  2  _  y+2 

rk+l  /rk+l  ri 


rk  +  /rk  -  ri 


Let  the  zone  width  (i.e.,  step  size)  =  Ar  and  let  i  =  0,  1,  2 
Then 

rk^  =  [(k  +  1  )Ar]2 
r2  =  (kAr)2 
r?  =  (iAr)2 

Thus 


N-1  (zone) 


n  -  n~V=y  =  SF  <Sk+l  '  Sk>  '°9 


(k  +  1)  +  [(k  +  1 )2  -  1*] 


1/2 


k  +  (k2  -  i2)l/2 


(3) 


4. 1.2. 2  Bimolecular  Mixing  Flow  Relations 

This  section  describes  how  density  and  mass  concentration  fields 
were  calculated  from  knowledge  of  the  index  of  refraction  field. 

For  a  multi -component  gas  case,  the  index  of  refraction  becomes, 
n  - '  *  £  Vi *  p  £  Vi 


where  K  =  Gladstone-Dale  constant 

C  =  mass  concentration  or  mass  fraction 
p  =  density 

subscript 

i  =  pure  species  i 
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V 


For  helium  and  air  (subscripts  H  and  A,  respecti vely) 


n  -  1 


and  because  by  definition 


-*»■ 


CA  +  CH  =  1 


n 


1+1 


(1) 


From  the  perfect  gas  law  the  density  ratio,  p/pm,  may  be  obtained  as 


P  _  p  M 
Pc  p  T  M 

00  '  CO  00 

where  T  =  temperature 
p  =  pressure 
M  =  molecular  weight 
superscript 

=  mean  quantity 

A  simplifying  assumption,  which  in  any  event  is  non-crucial  to  the  final 
results,  is  that  P  =  pro  and  T  =  T^,  i.e.,  the  pressure  and  temperature  of 
the  comparison  and  test  scenes  are  the  same.  Then,  p/pa  =  M/M 

From  Dalton's  Law  of  partial  pressures, 

p  =  T  £  piRi  *  M7=  PT| 


ft 
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where  R  =  gas  constant  per  unit  mass 
R  =  gas  constant  per  mole 
and  by  definition  M  is 


For  helium  and  air, 


Consequently  the  density  ratio  becomes 


The  case  here  is  for  a  comparison  scene  (sub  -)  comprised  of  air,  thus 


M»  =  mA  and  K„  =  ka 


and 


Eliminating  p/p^  fr0m  equations  (1)  and  (2)  one  obtains,  after  some  algebra 

("oo  ■  n) 


The  quantity  Kj>m  in  equation  (3)  is  linearly  dependent  on  pressure, 
and  only  slightly  dependent  on  temperature  and  the  wavelength  of  light 
used.  A  formula  established  by  Biot  and  Arago  expressing  these  dependencies 
and  confirmed  by  more  recent  experiments  (Reference  8)  is 
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nQ  -  1 

P«  K00  =  noo  ■  1  1  +  a  T  Poo 

00  00 

where  n  is  the  index  of  refraction  at  T  =  0°C,  p  =  1  atm 

is  the  coefficient  of  thermal  expansion  of  the  ambient  gas  in  C0_1 
is  the  ambient  temperature  in  C° 
p^  is  the  ambient  pressure  in  atmospheres 

Values  of  nQ  and  aK  required  for  this  study  were  obtained  from  tables  in 
Reference  8. 

The  above  analysis  pertains  to  an  air/helium  mixture.  When  N2  bubbles 
are  studied  in  an  SF6  environment,  the  subscripts  H  and  A  are  simply 
replaced  with  subscripts  N  and  S,  respectively,  and  subscript  °°  and 
subscript  S  are  identical.  Values  of  all  constants  used  are  tabulated 
below. 


Gases 

m 

Ma/Mh  or  Ms/Mn 

a 

00 

VKA  or  VKS 

8 

8 

CL 

i  Air/He 

1/2 

7.25 

0.893 

1.33  •  10"* 

1 

7.25 

1 

0.893 

2.67  ‘  10-4 

4 

7.25 

1.092 

0.893 

1.068  •  lO-3 

8 

7.25 

1.092 

0.893 

2.136  •  10-3 

10 

7.25 

1.092 

0.893 

2.67  •  10-3 

sf6/n2 

10 

5.22 

1.092 

1.98 

7.15  •  lO-3 

4.1.3  Fireball  Wake  Interferometry  Data  Reduction 

In  addition  to  mean  radial  helium  concentrations  through  the  core 
of  the  vortex,  radial  helium  concentrations  at  various  positions  through 
the  fireball  wake  were  extracted  from  the  interferometry  data.  A  free 
coordinate  system  and  normalisation  only  with  respect  to  the  wake  edges  were 
selected  for  averaging  the  wake  data,  conr.cnsurate  with  the  approach  used 
in  the  analysis  of  the  data  in  the  vortex  body.  The  same  procedure  applied 
to  both  infinite  fringe  and  finite  fringe  interferograms. 
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The  coordin  es  were  established  in  the  following  manner.  An  axis 
approximating  the  center  of  the  wake  was  determined  by  eye  and  drawn 
1  ngitudinally  through  the  core  of  the  vortex.  The  wake  coordinate  origin 
became  defined  by  a  line  drawn  normal  to  the  wake  axis  through  the  vortex 
coordinate  origin  (where  it  was  defined  in  Section  4.1.1.1).  All  cuts 
through  the  wake  were  parallel  to  this  line  (normal  to  the  wake  axis)  and 

their  locations  were  based  on  fractions  of  vortex  diameter  D  below  the 
wake  origin. 

Fringe  numbers  along  these  cuts  were  recorded  in  the  same  manner  as 
described  in  Section  4. 1.1. 2.  The  fringe  profiles  for  the  same  cuts 
were  averaged  over  n  interferograms  using  the  edges  and  center  as  normalizing 
parameters,  which  resulted  in  symmetric  averaged  fringe  profiles  when  n  wf.s 
sufficiently  large.  Since  the  axial  symmetry  condition  was  then  satisfied, 
the  previously  discussed  Abel  inversion  of  the  fringe  shift  equation  could 
be  applied  to  calculate  radial  mass  concentrations. 

4.2  SIZE/RISE  DATA 

Shadowgraph  movies  were  utilized  in  this  investigation  primarily  to 
study  gross  motions  of  the  rising  vortex  such  as  size  and  rise.  Only  a 
limited  amount  of  this  information  could  be  obtained  with  interferometry 
since  that  diagnostic  is  limited  to  obtaining  data  at  one  instant  of  time 
in  each  event.  After  acquisition  and  processing,  each  shadowgraph  film  was 
projected  on  a  large  sheet  of  paper  where  the  outline  of  the  vortex  was 
traced  frame  by  frame  until  the  vortex  disappeared  from  the  field  of  view. 
Since  only  gross  flow  features  appear  on  a  shadowgraph  frame  (change  of  ilium- 
i nation  proportional  to  the  second  derivative  of  density  rather  than  directly 
with  density  as  is  the  case  with  interferometry ) ,  distinct  features  such  as 
.addle  point  and  maxima  could  not  be  distinguished,  but  the  outline  of  the 
edge  could  accurately  be  identified.  While  viewing  the  shadowgraph  movies, 
it  was  possible  to  quite  accurately  establish  a  radial  axis  visually 

through  the  center  of  each  event.  The  averaging  of  this  size/rise  data  is 
outlined  below: 


event 


Referring  to  Figure  10,  one  may  define  for  the  jth  frame  of  the  ith 
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z.j  =  rise  height  from  top  of  bubble  tube 

D^.  =  diameter  of  center  of  vortex 

n-j j  =  inclination  of  diameter  to  horizontal 

x0- •  =  horizontal  displacement  of  center  of  vortex  from  center  of 
J  bubble  tube 

D0.  =  initial  bubble  diameter 

j  =  1  is  the  first  frame  after  the  bubble  is  broken 

The  following  physical  parameters  are  averaged  over  n  events  to 
yield  average  quantities  for  every  j:h  frame: 


mean  normalized  rise  height 


mean  normalized  diameter 


mean  rotation 


mean  normalized  translation  of  center  of 
vortex  from  bubble  tube  axis 


Again,  for  a  symmetric  event,  nj  and  should  approach  zero  as  a  check. 

In  addition,  standard  deviations,  a  are  computed  for  each  parameter  where 
a  for  some  parameter  P  is  defined  as 


CTP. 

J 


i  y  (p.. 

n  ^  ij 
i=l 


1/2 


where 


P. 

J 


-iy  p.. 

n  'i j 

i=l 


These  same  relationships  applied  equally  well  to  the  interferometry  data  and 
enabled  data  to  be  averaged  at  one  point  in  time. 
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5.  EXPERIMENTAL  RESULTS 


A  discussion  of  results  obtained  in  this  subscale  experimental 
study  are  presented  in  the  following  sections.  These  results  include 
vortex  size/rise  comparisons  from  both  interferometry  and  shadowgraph 
data,  qualitative  comparisons  of  the  internal  structure  between  laminar  and 
turbulent  vortices  based  on  particle  tracking  photographs,  selected  ma~.. 
concentration  profiles  from  interferometry  data,  and  a  composite  iso¬ 
concentration  picture  which  presents  the  spatial  distribution  of  low 
density  gas  within  an  averaged  turbulent  laboratory  event. 

5.1  VORTEX  SIZE  AND  RISE  RESULTS 

Figure  11  presents  vortex  size/rise  data  obtained  at  eight  different 
test  conditions  which  provide  over  a  factor  of  forty  (40)  variation  in 
Reynolds  number*.  The  vortex  rise  position,  z,  (given  by  the  center  of  the 
cloud  relative  to  the  center  of  the  initial  bubble)  is  presented  as  a  function 
of  the  vortex  diameter,  D.  The  results  presented  have  been  averaged 
over  each  data  set  at  a  given  test  condition  (see  Section  4.2  for  the  averaging 
technique  used  for  calculating  size/rise  information)  and  made  nondimensional 
using  a  characteristic  length  equal  to  the  initial  bubble  diameter,  Dq. 

The  shadowgraph  data  (shown  in  solid  symbols)  provide  results  up  to 
a  rise  position  of  about  5  bubble  diameters  since  delineation  of  the  vortex 
edge  beyond  a  few  bubble  diameters,  especially  in  laminar  cases,  was 
impossible.  Additional  data  points  are  provided  at  higher  rise  positions 
by  sets  of  averaged  interferometry  data  (shown  in  open  symbols).  The 
results  show  a  consistent  trend  w>'th  Reynolds  number,  namely,  the  higher 
Reynolds  number  data  display  a  greater  diameter  at  a  given  rise  position. 

Of  particular  significance  is  that  the  size/rise  data  fall  within  two 
distinct  bands;  one  corresponding  to  a  high  Reynolds  number  or  "turbulent" 
group  of  data  and  the  other  to  a  low  Reynolds  number  or  "laminar"  group. 


♦Reynolds  number  is  defined  as  Re  =  PooUyDg/y^  where  uj  is  the  terminal 
velocity  of  a  non-entraining  sphere,  D0  is  the  initial  bubble  diameter, 
and  and  y^  are  the  ambient  density  and  dynamic  viscosity,  respectively. 
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The  turbulent  data  band  includes  all  data  acquired  with  helium  in  air  at 
tank  pressures  of  8  and  10  atmospheres  and  with  nitrogen  in  sulfur 
hexaflouride  (SF6)  at  10  atmospheres  covering  a  range  of  Reynolds  numbers 
from  10,400  to  27,600.  The  laminar  data  band  encompasses  all  data  acquired 
with  helium  releases  in  air  at  1/2  and  1  atmosphere  pressure  and  for  initial 
bubble  diameters  ranging  from  0.9  to  2  inches  which  provides  Reynolds  numbers 
ranging  from  650  to  3000.  Additional  data  were  acquired  at  an  intermediate 
Reynolds  number  of  5200  using  a  helium  release  in  air  at  a  tank  pressure 
of  4  atmospheres.  These  data  tended  to  fall  between  the  laminar  and  turbulent 
bands  and  consequently  were  not  added  to  the  figure  (except  for  an  interferometry 
data  point)  to  improve  its  clarity.  These  data  are,  however,  tabulated  with 
the  other  size/rise  results  in  Table  I.  Mean  rotational  and  translational 
data  (n  and  xQ)  at  all  test  conditions  were  found  to  be  negligible,  as 
expected. 

Also  shown  in  Figure  11  are  bars  on  the  laminar  and  turbulent  data 
which  represent  the  standard  deviations  in  the  size  and  rise  data  at  that 
stage  of  development.  A  complete  summary  of  the  size/rise  standard  deviations 
for  both  the  shadowgraph  and  interferometry  data  are  presented  in  Table  II. 

These  results  are  given  for  each  set  of  averaged  interferometry  data  and 
for  the  averaged  shadowgraph  data  at  rise  heights  of  2,  3,  and  4  bubble 
diameters.  In  analyzing  these  variational  data,  it  was  found  that  there  were 
no  consistent  trends  with  Reynolds  number  or  stage  of  vortex  development. 

Instead,  all  data  showed  a  variation  from  the  mean  of  10  percent  ±  5  percent 
with  few  exceptions.  In  the  exceptional  cases,  (shadowgraph  data  at  8 
atmospheres  and  interferometry  data  at  1/2  and  1  atmosphere  pressure)  the 
standard  deviations  were  within  20  percent  of  the  mean.  It  should  be 
emphasized  that  the  standard  deviations  presented  are  for  size  and  rise 
dimensions  only,  and  in  general,  standard  deviations  of  internal  vortex 
characteristics  such  as  the  maximum  fringe  shift  from  the  interferometry 
data  tended  to  be  somewhat  larger. 

In  thoroughly  reviewing  the  interferometry  data,  it  was  found  that 
consistent  variations  in  the  configuration  of  the  vortices  (as  viewed  from 
the  horizontal)  were  present  as  a  function  of  Reynolds  number.  An 
additional  subjective*  dimension  beyond  those  described  in  Section  4  1  was 

♦Subjective  in  that  the  base  (bottom)  of  the  vortex  is  in  some  cases 
difficult  to  establish  because  of  the  presence  of  the  wake. 
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obtained  from  each  interferogram  to  provide  an  estimate  of  the  height 
(top  to  bottom  dimension)  of  the  vortex,  not  including  the  wake.  An 
average  vortex  diameter- to-height  ratio  was  then  obtained  for  four  individual 
sets  of  interferometry  data  to  provide  a  first  order  means  of  determining 
the  additional  entrainment  which  occurred  in  the  turbulent  events  relative 
to  the  laminar  than  would  be  indicated  by  only  their  respective  diameters. 
These  ratios  are  presented  in  the  table  below. 

Average  Diameter-to-Height  Ratios  of  Events  at  Different  Test  Conditions 

Reynol ds 
Number 

27,600 

13,000 

5,200 

1,300 

The  results  in  the  table  indicate  that  the  average  diamet^^-to-height 

ratio  for  the  turbulent  events  is  about  20  percent  less  than  for  the  laminar 

events  or  that  the  average  volume  of  the  turbulent  events  is  about  20  percent 

greater  than  the  volume  of  the  laminar  events  for  a  given  diameter.  This 

difference  in  configuration  results  in  a  greater  difference  in  entrainment 

between  the  laminar  and  turbulent  results  than  indicated  by  their  respective 

size/rise  characteristics  in  Figure  11.  In  terms  of  an  entrainment 
•  •  dr 

coefficient,  a  =  ^-  ,  the  rc  e  of  change  of  vortex  radius  with  rise  height 
(see  discussion  in  Section  2  of  Reference  1  for  example)  the  difference 
between  the  turbulent  and  laminar  data  would  be  increased  by  about  6  percent. 

A  surprising  result  is  the  much  greater  flattening,  or  larger  diameter- 
to-height  ratios  of  the  events  in  the  4  atmosphere  group  of  data.  This 
change  in  configuration  could  possibly  be  associated  with  the  transitional 
nature  of  the  flow  since  it  is  well  known  that  many  dramatic  effects  can 
occur  in  flow  fields  that  are  undergoing  transition  from  a  laminar  to  a 
turbulent  condition.  The  net  effect  is  that  on  a  volume  entrainment 
basis  the  4  atmosphere  set  of  data  tends  to  shift  and  collapse  with  the 
laminar  data 


Test  Condition 

SF6/N2,  10  Atmospheres,  Dq  =  1/2  Inch 
Air/He,  10  Atmospheres,  Dq  =  1  Inch 
Air/He,  4  Atmospheres,  Dq  =  1  Inch 
Air/He,  1  Atmosphere,  Dq  =  1  Inch 


Averaged  laboratory  fireball  rise  history  data  are  presented  in 

Figure  12.  The  same  nondimensional  rise  position,  Z/Dq  as  presented  in 

Figure  11  is  plotted  as  a  function  of  a  nondimensional  time  parameter, 

r  =  Vjj—  t  where: 
o 

t  is  the  fixed  time  for  the  events  to  reach  an  average  rise 
height  from  the  time  the  bubbles  were  broken 

g  is  the  gravitational  acceleration 

D  is  the  initial  bubble  diameter 
o 

Again  it  can  be  seen  that  low  Reynolds  number  or  laminar  data  fall  within 
one  group  and  the  turbulent  higher  Reynolds  number  data  fall  into  a  distinctly 
different  lower  rise  rate  pattern.  Curve  fits  to  both  the  laminar  and 
turbulent  data  are  shown  in  the  figure  to  better  illustrate  the  more  rapid 
rise  rate  of  the  laminar  data.  This  faster  rise  rate  exhibited  by  the 
laminar  vortices  results  from  a  less  pronounced  entrainment  of  air  in  these 
events  than  in  the  turbulent  events  which  is  consistent  with  the  size/rise 
and  configuration  data  presented  earlier. 

Again  the  intermediate  Reynolds  number  shadowgraph  data  acquired  at 
4  atmospheres  pressure  were  not  plotted  to  improve  the  clarity  of  the  figure. 
In  general,  these  data  fall  above  but  close  to  the  turbulent  data  curve 
fit.  In  addition,  three  turbulent  interferometry  data  points  (two  with 
SFb/N2  at  10  atmospheres  and  one  with  Air/He  at  3  atmospheres  pressure)  at 
high  rise  positions  were  not  shown  so  the  scale  of  the  plot  could  be 
expanded.  These  results  follow  the  same  trend  as  the  turbulent  data 
presented. 

Further  evidence  of  the  marked  differences  in  the  structure  and 
nature  of  the  laboratory  laminar  and  turbulent  vortices  beyond  that  already 
presented  in  terms  of  size,  rise,  and  shape  characteristics  is  displayed  in 
two  typical  particle  tracking  photographs  in  Figures  13  and  14.  These 
multiple  exposure  photographs  of  evolving  particulate  seeded  vortices  are 
presented  for  comparison  purposes  to  illustrate  dramatically  the  effect  of 
turbulence  on  toroid  structure.  Figure  13  represents  a  laminar  vortex  taken 
at  a  pressure  of  one  atmosphere  with  an  initial  bubble  diameter  of  one  inch 
corresponding  to  a  Reynolds  number  of  approximately  1300.  The  flash 
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Laboratory  Fireball  Rise  History 
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Figure  13.  Multiple  Exposure  Photograph  of  Particulate  Seeded 
Laminar  Fireball 
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rate  was  5.0  flashes  per  second.  The  important  feature  of  this  photograph 
is  the  distinct  toroid  "jelly-roll"  structure  apparent  throughout 
the  entire  rise. 

Figure  14  displays  a  turbulent  event  acquired  at  a  pressure  of  eight 
atmospheres  with  the  remaining  test  parameters  identical  to  the  one 
atmosphere  photograph.  The  Reynolds  number  for  this  particular  event  was 
about  10,400.  In  this  case,  turb'ilont  mixing  has  completely  overwhelmed  all 
internal  structure  evident  in  the  laminar  event.  It  is  apparent  from  these 
and  similar  vortex  structure  photographs  that  a  "jelly-roll"  structure  as 
such  is  strictly  associated  with  a  larinar  event  and  will  not  exist  in  a 
full  scale  detonation. 


5.2  COMPARISON  OF  TURBULENT  LABORATORY  DATA  WITH  FULL  SCALE  EVENT 

It  has  been  concluded  that  the  higher  Reynolds  number  (Re  >_  104) 
laboratory  data  are  indeed  turbulent.  Since  turbulent  flows  usually  exhibit 
a  weak  dependence  on  Reynolds  number,  one  would  hope  to  find  that  the 
laboratory  data  could  be  scaled  directly  to  full  scale  atmospheric  events. 

Only  two  characteristic  lengths  present  themselves;  1)  the  initial  diameter 
of  the  fireball  and  2)  the  atmospheric  scale  height.  Since  the  laboratory 
experiments  are  not  expected  to  be  sensitive  functions  of  the  latter,  one 
might  hope  that  the  initial  diameter,  Dq,  is  sufficient  for  scaling. 

Figure  15  shows  rise  versus  size  data  nondimensional  with  Dq.  The  laboratory 
data  which  have  factors  of  2  variation  in  the  initial  size  are  nicely  collapsed 
into  a  single  curve  as  are  the  data  of  the  full  scale  event  whose  initial 
diameter  is  approximately  10  times  as  large  as  the  laboratory  simulation. 

The  results  would  appear  to  confirm  that  Dq  is  the  only  important  characteristic 
length  and  that  the  Reynolds  number  variation  between  full  scale  and  laboratory 
conditions  is  not  critical  at  least  insofar  as  this  type  of  gross  dimensional 
data  is  concerned. 


For  the  time  dependence,  the  only  obvious  scaling  is  through  the 
gravitational  constant,  g,  which  when  coupled  with  the  initial  diameter, 

Figure  16  shows  D/Dq  plotted 

t  =V-^VjT-  t.  The  square  root  of  the  density 


Dq,  can  be  made  into  a  characteristic  time. 


against  a  nondimensional  time 
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Figure  15.  Comparison  of  Laboratory  and  Full  Scale  Fireball  Size/Rise  Data 


Figure  16.  Comparison  of  Laboratory  and  Full  Scale  Fireball  Size  History  Data 


ratio  has  been  included  to  account  for  the  density  variation  resulting  from 
the  altitude  change  that  the  full  scale  event  experiences  (Reference  9). 

This  ratio  is  identically  1  for  all  the  laboratory  data.  All  of  the  laboratory 
data  are  nicely  correlated  with  these  parameters  and  while  the  full  scale 
data  is  slightly  (^10%)  displaced,  nevertheless  the  trend  is  well  established. 
In  addition  to  agreeing  with  one  another,  all  of  the  data  are  well  represented 

2/5 

by  the  formula  D/Dq  ^  x  which  can  be  shown  to  be  a  result  of  turbulent 
entrainment  and  the  slowly  varying  toroidal  geometry  of  the  phenomena 
(Reference  9). 

It  is  concluded  that  the  laboratory  data,  when  viewed  in  the  proper 
nondimensional  frame  of  reference  are  representative  of  a  full  scale  atmos¬ 
pheric  event,  and  reasonable  credence  can  be  placed  in  the  concentration 
fields  thus  determined  which  are  discussed  in  the  subsequent  sectiois. 

5.3  LABORATORY  FIREBALL  MASS  CONCENTRATION  RESULTS 

Four  main  topics  are  discussed  in  this  subsection.  A  typical  averaged 
fringe  shift  distribution  is  described  which  satisfies  the  axisymmetry 
requirement  imposed  by  the  data  reduction  method.  Though  Reynolds  number 
dependence  on  mass  concentration  was  not  a  primary  objective  of  this  study, 
a  brief  comparison  between  laminar  and  turbulent  results  is  shown.  Next, 
mncentfation  distribution  dependence  on  rise  position  (or  equivalently, 
time  after  bubble  bursting)  and  the  geometric  similarity  that  results  are 
discussed  in  detail.  Finally,  mass  concentration  results  within  the  laboratory 
fireball  wake  are  discussed. 

5.3.1  Results  of  Data  Averaging 

In  Section  4.1.1,  the  technique  used  in  this  study  to  obtain  a 
composite  symmetric  event  by  ensemble  averaging  a  number  of  typically 
asymmetric  fringe  shift  distributions  was  described.  The  mean  distribution 
so  obtained  must  satisfy  the  inherent  requirement  present  in  the  data  reduction 
of  axisymmetry  before  it  may  be  used  to  compute  the  density  and  mass 
concentration  fields. 
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A  representative  mean  fringe  distribution  obtained  using  the 
aforementioned  averaging  technique  appears  in  Figure  17.  The  solid  curve 
corresponds  to  the  averaged  left  hand  side,  and  the  dashed  curve  to  the 
mean  right  hand  side.  The  particular  test  conditions  at  which  these  data 
were  acquired  were  as  follows:  tank  pressure  of  Id  atmospheres,  0.5  inch 
diameter  initial  N2  bubble  in  an  SF6  environment,  and  a  10.4  initial 
bubble  diameter  rise  1.175  seconds  after  bubble  breaking.  The  curves  result 
from  cuts  taken  through  the  center  of  the  vortices. 

It  is  seen  by  comparing  the  two  curves  that  the  averaged  event  is 
indeed  reasonably  axisymmetric.  Edge  y,adii  and  the  radial  locations  of  the 
fringe  maxima  deviate  by  less  than  2  percent.  The  relative  deviation  in 
the  fringe  shift,  S,  is  less  than  10  percent  at  most.  Prior  to  computing 
mass  concentration  profiles  through  an  averaged  event,  the  left  and  right 
hand  side  fringe  shift  distributions  were  averaged  together  to  provide  even 
a  b°tter  representation  of  a  symmetric  event. 

The  discrete  points  appearing  on  Figure  17  correspond  to  raw 
fringe  data  taken  directly  from  one  of  the  twenty-nine  interferograms  which 
compose  this  particular  set  of  data.  These  points  are  indicative  of  the 
asymmetries  present  in  indiviaiPl  events  and  the  necessary  utilization  of 
an  averaging  technique.  It  becomes  readily  evident  from  the  differences 
between  the  two  sides  that  an  Abel  inversion  of  fringe  data  from  single 
events  would  produce  misleading  results. 

5.3.2  Effect  of  Reynolds  Number  on  Averaged  Mass  Concentration  Profiles 

Figure  18  compares  a  laminar  mass  concentration  profile  with  one  that 
is  turbulent.  Both  profiles  were  reduced  from  data  acquired  at  a  common 
rise  position  and  have  their  radial  coordinates  nondimensional  ized  with  the 
initial  bubble  radii  for  comparison  purposes.  The  laminar  curve  evolves 
from  data  recorded  at  test  conditions  of  one  atmosphere  pressure,  a  one 
inch  initial  diameter  helium  bubble  in  an  air  environment,  and  at  0.55 
seconds  after  bubble  bursting.  The  turbulent  data  were  acquired  at  10 
atmospheres  pressure,  using  SF6/N2  as  the  gases,  a  1/2  inch  diameter 
initial  bubble,  and  a  time  after  bursting  of  1.175  seconds.  Both  data 
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Figure  17.  Averaged  and  a  Typical  Radial  Fringe  Shift  Distribution 


ar  and  Turbulenl 


sets  were  acquired  at  a  rise  position  of  approximately  6  initial  bubble 
diameters.  The  Reynolds  number  of  the  turbulent  data  (27,600)  is 
approximately  20  times  that  of  the  laminar  (1300). 

Two  main  features  become  evident  from  Figure  18.  First  is  that 
the  average  laminar  vortex  is  relatively  smaller  than  that  of  the  turbulent. 
Secondly,  the  stage  of  toroidal  development  of  the  turbulent  case  appears 
more  advanced  than  that  of  the  laminar.  This  is  evidenced  by  the  presence 
of  a  holp  at  the  center  of  the  turbulent  vortex  (see  Appendix  C).  The 

turbulent  toroid  also  exhibits  more  steep  concentration  gradients  near 
the  inner  edge. 

The  fact  that  the  turbulent  profile  has  a  lower  peak  mass  concentration 
than  the  laminar  is  inconclusive.  One  cannot  make  a  direct,  one-to-one 
comparison  between  magnitudes  since  the  initial  density  ratio  of  the 
laminar  data  (ratio  of  the  density  of  pure  helium  to  pure  air)  is  0.138 
whereas  that  of  the  turbulent  data  is  0.192.  These  initial  conditions,  if 
anything,  would  tend  to  reduce  the  differences  in  the  peak  magnitudes. 

For  the  sake  of  completeness,  this  comparison  of  laminar  and 
turbulent  mass  concentration  profiles  was  presented  and  discussed.  Strong 
evidence  that  the  turbulent  data  more  closely  resemble  the  full  scale  event, 
at  least  in  structure  and  gross  size/rise  character  (Section  5.2), 
motivates  placing  most  of  the  emphasis  in  the  remainder  of  Section  5  on 
the  turbulent  results,  and  only  when  additional  insight  may  be  gained  will 
laminar  results  also  be  presented.  To  be  complete,  all  data  acquired  during 
the  course  of  the  experiment  appears  in  tabular  form  in  Appendix  E. 

5’3-3  Hfect  of  Rise  Position  (or  Time)  on  Mass  Concentration 

A  more  conclusive  comparison  of  mass  concentration  data  than  that 
presented  in  Section  5.3.2  results  from  data  acquired  at  the  same  Reynolds 
number  but  at  different  rise  positions  (i.e.,  or  different  times).  Three 
large  repeat  sets  of  SF6/N2  10  atmosphere  data  were  acquired  for  this 
purpose  at  6.0,  10.4,  and  16.9  initial  bubble  diameters  rise,  or  0.625, 

1.175,  and  2.50  seconds  after  bubble  bursting,  respectively. 
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Figure  19  displays  central  cut  mass  concentration  results  for  the 
data  sets  mentioned  above.  The  radii  have  been  nondimensional ized  with 
the  initial  bubble  radius,  and  all  three  curves  have  been  plotted  to  the 
same  radial  scale  so  that  relative  vortex  growth  is  illustrated.  No 
attempt  was  made  to  1ocate  vertical  position  of  each  plot  to  the  scale 
of  the  horizontal  dimension. 

In  this  figure,  the  outer  radius  of  the  torus  grows  with  rise 
position  at  a  decreasing  rate  as  rise  position  or  time  increases.  This 
feature  was  evidenced  in  the  size/rise  discussions  and  need  not  be  amplified 
here.  More  importantly,  the  peak  concentrations  decrease  from  0.66  percent 
to  0.37  percent  to  0.21  percent,  their  peaks  become  flatter,  and  their 
inner  and  outer  concentration  gradients  decrease  with  increasing  rise. 

Data  pertaining  to  the  growth  of  the  inner  torus  radius  (i.e.,  the  hole  in 
the  center)  is  inconclusive,  because  what  defines  "hole"  is  somewhat 
ambiguous  since  the  data  reduction  accuracy  decreases  with  radius  (see 
Appendix  C). 

Figure  20  compares  the  same  data  which  appeared  in  Figure  19  but  in 
a  slightly  different  manner.  Each  mass  concentration  profile  has  been 
normalized  to  its  own  maximum  concentration,  and  the  radial  coordinate 
has  been  normalized  to  the  individual  edge  radius.  In  this  manner,  all 
three  curves  begin  with  a  zero  concentration  ratio,  increase  to  a  ratio  of 
one  at  its  maximum,  and  decrease  to  zero  at  rYrec|ge  =  1.  Figure  21  presents 
two  different  sets  of  laminar  data  in  an  identical  fashion.  These  particular 
laminar  data  were  acquired  at  one  atmosphere  pressure,  with  air/helium  as 
gases,  and  at  rise  positions  of  3.1  and  5.9  initial  bubble  diameters,  res¬ 
pectively. 

A  comparison  of  the  three  curves  in  Figure  20  shows  that  there  is 
very  little  difference  at  all  between  them,  which  is  indeed  a  most  interesting 
finding.  This  indicates  that  for  all  practical  purposes,  the  mass 
concentration  distribution  in  the  torus  remains  geometrically  similar  after 
torus  formation  has  been  established  for  this  range  of  rise  positions,  in 
the  laminar  case  (Figure  21),  though  only  two  sets  of  data  were  analyzed  and 
at  comparatively  lower  rise  positions,  the  same  results  bear  credence. 

4  V 
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gure  19.  Averaged  Radial  Mass  Concentration  at  Various 
Positions;  Turbulent  Results 


Figure  20.  Normalized  Mass  Concentration  Distributions 
Turbulent  Data  at  Various  Rise  Positions 


5.3.4  Laboratory  Fireball  Wake  Results 

Two  sets  of  data,  one  turbulent  and  one  laminar,  were  analyzed  to 
investigate  the  character  of  the  laboratory  fireball  wake  and  to  provide 
insight  as  to  the  significance  of  the  wake  in  relation  to  the  main  vortex 
flow  field. 

An  averaged,  10  atmosphere  SF6/N2  vortex  which  rose  10.4  initial 
bubble  diameters  provided  the  data  base  for  the  turbulent  wake  calculations. 
Three  cuts  were  made  through  the  wake  in  the  fashion  outlined  in  Section 
4.1.3  at  3/10  D,  1/2  D,  and  3/4  D  below  the  defined  fireball  coordinate 
origin.  A  physical  constraint,  namely  the  format  size  of  the  interferogram, 
precluded  any  analysis  of  data  below  3/4  D.  The  laminar  wake  results  were 
ascertained  from  an  averaged,  one  atmosphere,  air/He  set  of  interferograms . 
Since  the  laminar  vortices  were  smaller  than  the  turbulent  ones,  results 
were  calculated  along  cuts  at  1/2  D,  3/4  D,  and  D  normal  to  the  "wake 
axis. " 

Figure  22  provides  nondimensional  comparison  between  the  turbulent 
results  and  those  of  the  laminar.  Mass  concentrations  along  the  various 
cuts  through  the  wake,  Cw>  are  divided  by  the  maximum  concentration  along  the 
central  cut  of  its  corresponding  average  vortex,  Cv  .  The  radial  position 

Mia  A 

along  the  wake  cut,  r  ,  is  nondimensional ized  by  the  initial  bubble  radius  r 

W  0 

to  suppress  the  effects  of  different  initial  bubble  sizes  (0.5  inch  diameter 
for  the  turbulent,  0.4  inch  for  the  laminar). 

There  are  several  wake  characteristics  that  are  illustrated  in 
Figure  22  worthy  of  note.  First,  the  turbulent  wake,  both  relative  to  the 
initial  bubble  size  and  to  its  corresponding  vortex  size,  is 
considerably  wider  than  the  laminar  wake.  Secondly,  the  peak  wake  concentra¬ 
tion  of  the  turbulent  case  is  far  less  than  that  of  the  laminar,  and  remains 
reasonably  constant  throughout  much  of  the  width  of  the  wake  whereas  the 
concentration  of  the  low  density  gas  in  the  laminar  wake  decreases  rather 
steeply  from  its  maximum  to  the  edge. 
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A  calculation  presented  in  Appendix  D  and  discussed  more  thoroughly 
in  the  following  section  indicates  that  slightly  less  than  10  percent 
of  the  nitrogen  present  in  the  vortex  flow  field  is  contained  in  the 
wake.  To  the  degree  of  accuracy  of  the  calculation  made,  this  indicates  that 
most  of  the  initial  mass  does  indeed  remain  in  the  vortex  after  torus 
formation  has  been  established.  One  could  reasonably  extrapolate  this 
result  to  a  full  scale  nuclear  event  and  conclude  that  the  assumption 
that  all  weapon  and  carrier  debris  remains  in  the  torus  would  appear  to 
be,  in  the  first  approximation,  valid. 

5.4  COMPOSITE  ISO-CONCENTRATION  PICTURE  FOR  A  TURBULENT  VORTEX 

Iso-concentration  profiles  for  a  turbulent  vortex  were  derived  from 
a  series  of  averaged  concentration  profiles  and  are  presented  in  Figure  23. 
The  set  of  interferometry  data  which  was  used  in  preparation  of  this  figure 
was  obtained  with  an  N2  release  in  SF6  at  10  atmospheres  pressure  after  the 
center  of  the  events  rose  a  distance  of  10.4  initial  bubble  diameters.  The 
individual  concentration  profiles  through  the  vortex  body  and  wake  are 
presented  in  Appendix  E  since  they  may  be  more  useful  in  that  form  for 
comparison  to  code  results.  The  concentration  level  corresponding  to 
each  iso-concentration  profile  is  shown  in  Figure  23.  The  grid  marks 
shown  along  the  axes  are  0.25  inches  apart  which  corresponds  to  the  initial 
bubble  radius  of  this  averaged  event. 

A  calculation  of  the  mass  of  nitrogen  in  the  portion  of  the  event 
shown  was  made  (see  Appendix  D  for  details)  which  yielded  a  result  equal  to 
about  87  percent  of  the  original  N2  in  the  unbroken  soap  bubble.  It  should 
not  be  inferred  that  13  percent  of  the  original  mass  is  necessarily  in  the 
wake  below  the  averaged  event  shown  since  the  number  of  profiles  used  and 
the  averaging  necessary  in  the  integration  technique  was  only  accurate  to 
within  about  10  percent.  It  was  found  that  slightly  over  90  percent  of  the 
calculated  mass  of  N2  was  located  in  the  vortex  body  or  less  than  10  percent 
of  the  total  was  located  in  the  wake  shown  in  Figure  23. 
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Iso-c ncentration  Profiles  for  a  Turbulent  Vortex 
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In  rising  over  10  bubble  diameters  the  averaged  vortex  shown  has 
entrained  a  mass  of  SF6  more  than  600  times  greater  than  the  mass  of  N2 
(Appendix  D) .  As  a  result  of  this  entrainment  the  iso-concentration  lines 
range  from  a  peak  of  slightly  over  0.35  percent  at  the  center  of  the 
vortex  core  to  zero  at  the  edge.  It  should  be  noted  that  the  edge  profile 
shown  as  a  zero  concentration  line  is  an  estimate  only  since  defining  the 
exact  location  of  zero  fringe  shift  is  impossible.  It  is  also  interesting 
to  note  that  the  region  defined  by  concentrations  above  0.1  percent  consists 
of  a  toroid  with  a  kidney-shaped  cross-section. 

The  largest  mean  concentration  gradients  can  be  seen  to  be  on  the 
outside  top  nd  sides  of  the  vortex  and  on  the  inside  of  the  torus  near  the 
core.  Again  it  should  be  noted  that  the  mean  concentration  profiles 
presented  here  have  resulted  from  the  averaging  of  29  events  and  consequently 
instantaneous  gradients  at  the  edge  of  a  single  turbulent  event  would 
undoubtedly  be  greater. 

A  transparent  overlay  shows  both  a  sphere  which  encompasses  the 
bulk  of  the  averaged  event  and  a  torus  which  provides  a  reasonable  fit  to 
the  higher  concentration  level  torus  core.  It  can  be  seen  that  a  sphere 
does  provide  a  rough  but  fairly  good  representation  of  the  volume 
contained  in  the  vortex  flow  field  with  the  primary  exceptions  being  the 
depression  and  hole  in  the  top  and  the  wake  structure  aft  of  the  torus. 

The  torus  shown  in  the  overlay  is  geometrically  similar  to  that  described 
by  geometric  d-Hca  from  the  same  low  yield-low  altitude  detonation 
discussed  earlier  in  connection  with  Figures  15  and  16. 


6.  CONCLUSIONS 


Listed  below  are  the  conclusions  considered  significant  based  upoi 
results  presented  in  Section  5. 

1.  Mean  size/rise  characteristics  of  turbulent  and  laminar 
events  are  distinctly  different  and  fall  into  two  well- 
defined  bands.  The  turbulent  data  display  a  greater  size 
for  a  given  rise  distance. 

2.  The  larger  diameter  and  different  shape  of  the  turbulent 
events  reflect  an  increased  rate  of  entrainment  over  that 
which  exists  in  the  laminar  cases. 

3.  As  a  result  of  a  decreased  rate  of  entrainment,  the  laminar 
data  exhibit  a  greater  rise  rate  than  the  turbulent  data. 

4.  Particle  tracking  photographs  show  a  "jelly-roll"  structure 
in  toroids  at  low  Reynolds  number  laminar  conditions. 

Turbulent  mixing  and  large  scale  turbulent  structure  completely 
overwhelm  any  wrap-up  structure  at  the  higher  Reynolds 

number  test  conditions. 

5.  Laboratory  turbulent  size/rise  data  are  in  excellent  agreement 
with  data  from  a  full  scale  detonation  when  results  are  viewed 
in  the  appropriate  nondimensional  frame  of  reference. 

6.  A  nondimensional  time  parameter  (x  =  %—  Va=t)  has  been 

P  Do 

shown  to  correlate  laboratory  and  full  scale  fireball 
growth  rate  data.  The^e  data  show  that  the  fireball  diameter 
is  proportional  to  . 

7.  Comparison  of  events  at  an  identical  Reynolds  number  but  at 
different  rise  positions  indicates  that  the  peak  concentration 
levels  decrease,  the  peaks  become  less  pronounced,  and  the 
concentration  gradients  decrease  with  increasing  rise  distance. 

8.  Concentration  distributions  in  a  rising  vortex  have  for  all 
practical  purposes  been  shown  to  be  geometrically  similar  after 
the  toroid  has  become  established  (for  a  range  of  rise 
positions  from  6  to  17  initial  diameters). 

9.  A  distinct  wake  is  present  below  a  rising  toroid.  However,  it 
was  found  that  at  a  high  Reynolds  number  condition,  less  than 

10  percent  of  the  total  mass  was  in  the  wake,  or  over  90  percent 
was  present  in  the  vortex  body. 
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10.  At  the  same  high  Reynolds  number  condition,  it  was  found  that 
a  vortex  entrained  more  than  600  times  its  original  mass  in 
rising  a  distance  of  about  10  initial  bubble  diameters. 

11.  The  largest  mean  concentration  gradients  in  an  averaged 
turbulent  vortex  were  found  to  be  on  the  outside  top  and  sides 
of  the  vortex  and  on  the  inside  of  the  toroid  near  the  core. 

12.  It  has  been  found  that  for  the  same  high  Reynolds  number 
condition  that  the  higher  concentration  levels  were  located 
in  a  small  kidney-shaped  core. 


7.  RECOMMENDATIONS 


The  laboratory  fireball  data  described  in  this  report  have  provided 
for  the  first  time  a  quantitative  data  base  in  terms  of  mean  mass  concentra¬ 
tion  and  mean  density  field  distributions  for  use  in  the  checkout  of 
hydrodynamic  codes.  Mean  shadowgraph  data  along  with  interferometry  data 
were  ensemble  averaged  to  provide  laboratory  size/rise  characteristics. 

The  high  Reynolds  number  turbulent  laboratory  size/rise  data  are  in  very 
good  agreement  with  data  from  a  low  yield,  low  altitude  nuclear  event  which 
demonstrates  that  this  experimental  technique  does  simulate  the  important 
physics  that  controls  fireball  hydrodynamics.  Strong  evidence  that  the 
simulation  is  credible  coupled  with  the  need  for  quantitative  data  to 
support  code  calculations  indicate  that  further  subscale  laboratory  studies 
are  justified.  Data  similar  to  that  presented  in  this  report  could  be 
produced  in  a  very  straightforward  manner  using  the  existing  experimental 
facility  and  data  analysis  techniques.  Direct  turbulence  measurements 
could  be  made  through  the  addition  of  two  primary  diagnostic  instruments. 

Additional  mean  concentration/density  data  could  be  obtained  without 
facility  modification  or  development  of  new  data  analysis  techniques.  If 
further  support  is  required  for  comparison  to  code  results,  more 
iso-concentration  pictures  (e.g.,  Figure  23)  could  be  generated  for 
various  stages  of  vortex  development.  In  addition,  these  pictures  would 
provide  better  physical  insight  into  the  flow.  These  unique  data, 
quantitative  in  nature,  could  be  generated  at  either  laminar  or  turbulent 
conditions,  depending  on  the  code  development  requirements. 

An  extension  of  the  study  would  provide  quantitative  turbulence 
measurements,  necessary  to  guide  and  provide  inputs  for  the  development  of 
turbulent  model  equations  for  fireball  flow  fields.  Direct  density 
fluctuation  and  velocity  fluctuation  measurements  could  be  made 
simultaneously  with  the  addition  of  two  diagnostic  instruments.  The  density 
fluctuation  measurements  could  be  made  using  a  small,  fast  response  probe 
developed  and  successfully  used  at  California  Institute  of  Technology  by 
Brown  and  Rebollo  (Reference  12).  A  dual  scatter,  laser  Doppler  velocimeter 
(LDV)  as  used  by  Sullivan,  et.  al  (Reference  13)  could  be  incorporated  into 
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the  facility  in  order  to  measure  two  orthogonal  velocity  components  of  the 
vortices  simultaneously.  With  these  measured  fluctuation  quantities,  the 
important  turbulent  transport  terms  in  the  Reynolds  averaged  momentum  and 
thermal  energy  (and/or  species)  conservation  equations,  namely  the  turbulent 
Reynolds  stresses  and  the  turbulent  heat  or  species  fluxes  are  known.  Thus, 
a  two-dimensional  transient  solution  of  these  conservation  equations  may  be 
accomplished. 


r 
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Table  I 


Mean  Vortex  Size/Rise  Data 


Shadowgraph  1  Atm. 

2"  Bubble  Re  =  3000 


t  (sec) 

Z/°o 

d/Dq 

X/Dn 

0  0 

n  (deg) 

.0417 

.064 

1.130 

-.006 

.06 

.0834 

.600 

1.223 

.002 

1.79 

.1251 

1.044 

1.369 

.041 

2.75 

.1668 

1.422 

1.525 

.056 

.23 

.2085 

1.782 

1.654 

.059 

1.52 

.2502 

2.049 

1.712 

.050 

3.29 

Shadowgraph  4  Atm. 

1"  Bubble  Re  =  5200 


t  (sec) 

Z/Do 

D/Dq 

VDo  r 

i  (deg) 

.0417 

.251 

1  . 

,385 

-.008 

-1.35 

.0834 

.632 

1  , 

.540 

-.020 

-1.00 

.1251 

1.037 

1. 

,606 

-.028 

-1.37 

.1668 

1.440 

1  , 

,729 

-.026 

-.93 

.2085 

1.804 

1  . 

,831 

-.030 

-1.19 

.2502 

2.217 

1. 

,958 

-.020 

-.20 

.2919 

2.543 

2. 

,076 

-.036 

-.04 

.3336 

2.915 

2. 

,225 

-.011 

.51 

.3753 

3.278 

2, 

,275 

-.031 

.64 

.417 

3.593 

2 

.312 

-.137 

.14 

Shadowgraph 

8  Atm. 

1 

1 

Shadowgraph 

10  Atm. 

1"  Bubble 

Re  =  10 

o 

o 

vr 

1"  Bubble 

Re  =  13,000 

t  (sec) 

z/Do 

D/Do 

V°o 

n  (deg) 

t  (sec) 

z/Do 

D/Do 

VDo  * 

(deg) 

.04V 

.138 

1.418 

-.001 

.20 

.0417 

.286 

1.301 

-.002 

-.35 

.0834 

.554 

1.622 

-.004 

.37 

.0834 

.686 

1.489 

.008 

-.08 

.1251 

1.038 

1.830 

-.008 

.38 

.1251 

1.178 

1.717 

.04 

.47 

.1668 

1.513 

2.003 

.009 

.83 

.1668 

1.598 

1.974 

.018 

.66 

.2085 

1.927 

2.219 

-.003 

-.09 

.2085 

1.967 

2.128 

-.006 

.73 

.2502 

2.294 

2.316 

-.061 

.34 

.2502 

2.358 

2.305 

.005 

2.58 

.2919 

2.683 

2.434 

-.087 

1.41 

.2919 

2.731 

2.483 

.008 

4.71 

.3336 

3.022 

2.507 

-.077 

2.11 

.3336 

3.074 

2.577 

.030 

5.16 

.3753 

3.347 

2.575 

-.067 

2.12 

.3753 

3.441 

2.639 

.022 

5.06 

.417 

3.732 

2.722 

-.051 

1.84 

.417 

3.794 

2.774 

.054 

4.98 

.4587 

4.004 

2.813 

-.047 

2.92 

.4587 

4.142 

2.914 

.086 

4.67 

.5 

4.182 

3.026 

0. 

1.79 

.5 

4.426 

3.001 

.082 

4.20 

.5417 

4.490 

3.181 

.008 

2.2 

.5417 

4.505 

2.959 

.020 

5.10 

.5834 

4.600 

3.304 

.065 

3.11 

.5834 

4.730 

3.162 

.087 

4.72 

Table  I  (CTD) 

Mean  Vortex  Size/Rise  Data 


Interferometry  Data 


Gases 

P  (Atm.) 

Dq  (in.) 

Re 

t  (sec) 

z/Do 

D/D0 

Air/He 

1/2 

1.0 

650 

.48 

5.01 

2.88 

Air/He 

1 

1100 

.19 

3.16 

2.22 

Air/He 

1 

1.0 

1300 

.55 

5.92 

2.85 

Ai r/He 

4 

1.0 

5200 

.80 

5.65 

3.12 

Air/He 

8 

1 .0 

10,400 

1.44 

9.44 

4.84 

Air/He 

10 

1 .0 

13,000 

.90 

5.76 

4.04 

SVN2 

j  10 

0.5 

27,600 

.63 

5.99 

3.80 

sf6/n2 

10 

0.5 

27,600 

1.18 

10.39 

5.12 

SF6/N2 

10 

0.5 

27,600 

2.50 

16.8 

7.26 
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Table  II 

Vortex  Size/Rise  Data  -  Standard  Deviations  a 

Shadowgraph  Data  -  Various  Rise  Positions 
Z/D0  *  2 


Gases 

P  (Atm.) 

°0  (1n- 

)  Re 

°z 

aD/D 

Air/He 

1 

2.0 

3000 

2.049 

.149 

.120 

Ai r/He 

4 

1.0 

5200 

2.011 

.121 

.144 

Ai r/He 

8 

1.0 

10,400 

1.927 

.200 

.120 

Ai r/He 

10 

1  .0 

13,000 

1.967 

.139 

.120 

z/D0  *  3 

Air/He 

4 

1.0 

5200 

2.915 

.131 

.121 

Air/He 

8 

1.0 

10,400 

3.022 

.187 

.087 

Ai r/He 

10 

1.0 

13,000 

3.074 

.107 

.130 

z/DQ  %  4 

Air/He 

8 

1.0 

10,400 

4.004 

.170 

.124 

Ai r/He 

10 

1.0 

13,000 

4.142 

.083 

.  117 

Interferometry  Data 

Air/He  1/2 

1.0 

650 

5.01 

.137 

.177 

Ai r/He 

1 

0.9 

1100 

3.16 

.192 

.181 

Ai r/He 

1 

1.0 

1300 

5.92 

.028 

.128 

Ai r/He 

4 

1.0 

5200 

5.65 

.097 

.119 

Ai r/He 

8 

1.0 

10,400 

9.44 

.104 

.066 

Ai r/He 

10 

1.0 

13,000 

5.76 

.083 

.101 

sf6/n2 

10 

0.5 

27,600 

5.99 

.093 

.118 

SF6/N2 

10 

0.5 

27,600 

10.39 

.213 

.149 
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APPENDIX  A 

Characterization  of  Particle  Path  Dependence 


To  be  an  effective  fluid  tracer,  particles  must  possess  relatively  low 
inertias  and  high  drags  to  follow  the  streamlines  accurately.  At  the 
same  time,  thry  must  be  sufficiently  large  to  become  visiole  under  intense 
illumination.  Particle  tracking  is  particularly  valuable  to  describe  flow 
patterns  not  easily  accessible  to  other,  more  standard  methods  of  measurement. 
The  utility  of  this  method  depends  on  the  accuracy  of  streamline  tracing 
desired. 

Assuming  negligible  radial  velocities  and  net  external  forces,  and 
using  Oseen's  approximation  for  CD(Re) ,  results  from  an  analysis  by  Wright 
(Reference  10)  show  that  the  path  dependence  of  a  particle  in  a  vortex  may 
be  characterized  by  a  single  parameter, 


,  _  O  /OTrayD1  \ 

o  ~  V.  '  M'  } 
o 


where  Rq  =  initial  position  of  particle  from  center  of  core  of  vortex 
initial  fluid  speed 
radius  of  particle 
w  =  fluiu  viscosity 


o 

a 


D' 


M'  = 


PP  = 
P  = 


Cn^e 

drag  parameter  =  =  1  +  0.1875(Re)  -  0.0148( Re) 2  +  ... 

(Reference  11) 

appirent  mass  of  particle  =  j  Tra3(pp  +  -  p) 

density  of  particle 
fluid  density 


A  plot  of  this  particle  path  dependence  on  parameter  from 

Reference  10  is  presented  in  Figure  A.l.  It  is  evident  that  large  X 
is  desirable. 
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Dependence  on  Parameter 


For  the  particles  used  in  this  study,  a  value  of  X  may  easily  Lie 
computed.  An  example  is  presented  to  illustrate  what  happens  to  a  1  on  micron 
diameter  particle  within  a  vortex  having  a  core  radius  of  1  inch,  a 
circumferential  velocity  of  1  ft/sec,  and  the  fluid  predominantly  air  at 
1  atmosphere.  The  particle  Reynolds  number  for  this  case  is  calculated  to 
be  2,  thus  the  drag  parameter  O'  using  Oseen's  approximation  is  1.33.  For 
this  particular  case,  XQ  -  14  which  based  upon  Figure  A-l  is  indicative 
of  a  very  good  tracer.  Since  D'  increases  strongly  with  Re  whereas  M'  is 

very  weakly  dependent  on  Re,  the  particulate  becomes  a  better  tracer  in  the 
higher  pressure  regimes. 
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APPENDIX  B 


ABEL  INVERSION  OF  THE  AX I SYMMETRIC  FRINGE  SHIFT  EQUATION 
The  inversion  of  the  integral  equation 

s(,1,4r°<n‘nJrfi ™  (B'" 

for  change  in  index  of  refraction  (n  -  nj  expressed  in  terms  of  fringe  shift 
S  is  by  no  means  trivial,  and  the  details  of  this  inversion  form  the  scope 
of  this  appendix. 

Introduce  v  =  r2 

u  =  y2 

Then , 


(B-2) 


Multiply  both  sides  of  (B-2) 


d  u 

by  ZH  and  integrate  from  w  to  v  . 
Ai  -  w  0 


s(u)du  „  i  r\  du  r o  (n  -  nJdv 

/u  -  w  *  w  /u  -  w  u  /v  -  u 


03-3) 


Interchange  the  order  of  integration  on  the  right  hand  side  of  (B-3) 
in  the  following  manner  (note  the  change  in  limits) 


r°  ^-u-  =1  /‘V°  (n  -  n  )dv  f"0  dLI 
‘'w  /u  -  W  “W  °°  *W  /(u  -  w)  (v  -  u) 


(B-4) 
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v 

The  integral  / 

J  W 


w  /|U  ■  w)  (V  -  u) 
variable  of  integration  p  such  that 


may  be  evaluated  by  introducing  a  new 


u  =  v  +  p(w  -  v) 


Then  the  integral  reduces  to 


rvo _ du _ r 

w  -  w )(v  -  u)  JQ 
Consider  the  definition  of  the  Beta  Function 

B(m,n)  =  J  u"1”1  (i  -  u)n_1du 
o 

In  this  case,  from  (B-5)  it  is  seen  that  m  =  n 

r(l)  =  1, 


r(m)r(nj 
r(m  +  n] 


(B-5) 


\  ,  and  since  r(^-) 


=  /iF  and 


;u  -  w)  (v  -  u] 


Thus, 


W  •'1,1  T~_ 


'  u  -  w 


(B-6) 


Now,  n  -  n^  may  be  obtained  by  differentiating  bo.fch  sides  with  respect  to 
w  and  then  applying  the  Leibnitz  Rule. 

The  left  hand  side  of  (B-7)  becomes 


l  ("  -  nJdv  *  -  »jw 


(B-8) 
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Thp  right  hand  integral  may  be  integrated  by  parts 


f  °  -^-du  =  2  S(u)  AT 
J  w  /u  -  w 


w 


w 


°-2/s(v“v' 

yS(w) 


w  dS 


By  definition,  S(v  )  is  identically  0  so  the  first  term  is  0,  and 


letting  dS  =  —  du 


(n  ■  n->w  =  f  ®  /  V°  £  <*u  (8-9) 

To  evaluate  the  derivative  of  the  integral,  Leibnitz'  Rule  is  again  applied. 


d_ 

dw 


W  W 


/u  -  w  -r-  du 
du 


Therefore , 


dS 

du 


du 


(n  -  n  ) 

*  W  71  ^  /iT^ 


which  may  be  equivalently  written  as 


(n  -  n  )  —  f  0 

°°  r=y  tt  J 


r  dS  Hr 

r-  a?  dr 


(r2  -  y2)>/2 


which  is  Equation  (2)  in  Section  4.1.2. 


APPENDIX  C 


EFFECT  OF  ERROR  IN  ^  NEAR  THE  SINGULARITY  IN 
THE  ABEL  INVERSION  OF  THE  FRINGE  SHIFT  EQUATION 

Motivated  by  the  occasional  appearance  of  slightly  negative  mass 
concentrations  in  the  region  near  r  =  0  (these  were  omitted  in  the  plots 
but  included  in  the  tabulated  results  in  Appendix  E),  an  effort  was 
undertaken  to  ascertain  the  origin  of  these  physically  unrealistic  results. 
Consider  Equation  (3)  in  Section  4. 1.2. 2  with  subscripts  modified  for 
gases  SF6/N2  rather  than  air  and  helium. 


-(n  -  n  ) 

oo 


t 

P  K 

00  00 

N  /n  -  n  ' 

(  00  +  i 

UMs 

(Ku 

\p  K  1 

'  00  00 

/UN  -  v 

'W *  ) 

From  the  table  at  the  end  of  Section  4. 1.2. 2, 

,  MS  KN 

Poo  =  7.15  •  10-3;  ^  =  5.22;  f- =  1.98 

N  S 

For  the  10  atmosphere  10.4  bubble  diameter  rise  data,  it  was  found  that 

8  •  10-6  <  vn  -  n)  <  8  •  10-5 

Thus,  1 0” 3  <  ~ — —  <  10"2 

Poo  \o 

n  “  noo 

And  the  - — ^ —  term  in  the  denominator  becomes  negligible  compared  to  1. 

^00  oo 

Thus,  for  this  case 

CH  'v  (n  -  nj 
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Consider  the  equation 


n  n  \  -  zi  /  e  dr 

“  r=Y  *  Jy  (r2  .  y2)1/2 

Typical  values  of  ^  (which  in  the  data  reduction  are  approximated  by  — ■) 
in  the  region  near  the  singularity  are  *oproximately  half  a  fringe  in 
0.1  inches;  or  about  5  fringes/inch.  Suppose  that  locally  there  may  be 
an  error  of  ±20  percent  from  the  true  ^  ,  which  may  be  quite  likely  when 
one  measures  slope  off  of  the  reconstructed  interferograms. 

An  error  of  20  percent  of  5  fringes  per  inch  (i.e.,  1  fringe/ inch) 
would  modify  the  original  calculation  the  following  way: 


^  (&*  l)dr 

”•4  (r2-y2)l/2 


2k_ 

■n 


f 


L  y 


(r2  -  y2) 


T W 


dr 


(r2  .  yi)'!2 


Original  Calculation 


Error 


Once  the  integral  on  the  right  is  evaluated,  the  incurted  error  in  (n  -  nj 
becomes 


+  -  log 

TT 


This  indicates  that  the  error  in  Cu  introduced  by  an  error  in  4  of  1 

H  / 2  r  \  dr 

fringe/inch  is  directly  proportional  to  log  (-y— j  •  A  table  of  this  error  as 


a  function  of  r/rg  appears  below: 

/2r  \ 

y/re 

-5 

o 

l0»(-y£) 

0 

0  (center) 

00 

0.05 

.256 

3.689 

0.1 

.512 

2.996 

0.2 

1 .024  ('v  inner  edge) 

2.303 

0.5 

2.56  (<v  peak  CM) 

1.3  86 

1 .0 

5.12  (outer  edge) 

.693 
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For  very  small  y,  that  is,  for  a  calculation  made  for  a  point  near  the 
origin,  the  error  induced  by  an  error  in  slope  becomes  large  indeed. 

When  the  calculation  is  made  at  points  further  away  from  the  origin  (larger 
y),  the  error  is  less  pronounced.  It  was  found  that  an  error  in  slope  of 
±1  fringe  was  indeed  sufficient  to  make  the  calculated  negative  mass 
concentration  positive.  In  this  report,  negative  concentrations  near  the 
vortex  origin  are  ignored,  concentrations  in  that  region  are  just  considered 
negligible. 
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APPENDIX  D 


CALCULATION  OF  MASS  OF  A  LABORATORY  FIREBALL 
D-l .  MASS  OF  LOW  DENSITY  GAS 

This  calculation  was  undertaken  to  compare  the  amount  of  nitrogen 
present  in  an  analyzed  toroid  to  that  in  the  unbroken  bubble.  It 
accomplishes  two  purposes;  first,  it  tests  the  validity  of  Figure  23,  and 
second,  gives  the  relative  amounts  of  mass  in  both  the  vortex  and  the  wake. 

From  Dalton's  Law  of  Partial  Pressures, 


PN  PS 


p 


00 


where  PN  and  p$  are  the  partial  pressures  of  N2  and  SF6  respectively,  and 
pn  and  ps  are  the  densities  of  N2  and  SF6  in  the  disturbed  region.  This 
equation  utilizes  the  facts  that  pressure  and  temperature  are  constant 
everywhere  and  that  the  ambient  mediun  is  composed  of  pure  SFb.  Multiplying 
through  by  Mr,  we  get 


The  density  in  the  mixed  gas  region  is  equal  to  the  sum  of  the  densities 
of  N2  and  SF6  since  the  vol  tme  is  fixed.  Thus, 

pN  +  pS  =  p  (D-2) 

Replacing  P$  in  (D-l)  with  p  -  pN  fron  (D-2)  and  integrating  o\er  the  known 
volune,  the  mass  of  nitrogen  within  the  volume  becomes 


m 


N 


«s\  /  0  '  t)dV 


( D-3) 
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The  mass  of  the  spherical,  initially  pure  nitrogen  bubble  is 


X  =  3  »  r0  pN0 

where  rQ  is  the  bubble  radius  and  pN  is  the  density  of  pure  nitrogen. 
This  may  be  modified  to 


X  =  J  *  ro  ^  P. 


(D-4) 


The  ratio  between  the  calculated  mast  and  the  initial  mass  becomes 


1  '  4  "ro3 


/  0  -  MdV 


(D-5) 


Consider  now  the  inteqral  I  (1  -  S— iwu 

y  J  Poo  dV'  Break  UP  the  volume  into  M  thin 

<  "leal  volumes  each  having  thickness  Ay. .  Bearing  in  mind  that  the 

volume  of  any  cylinder  may  be  expressed  as  2„ay  f  *  rdr,  where  r  is  the 

edge  radius,  the  integral  in  question  may  be  expressed  in  terms  o/a  sum  of 
integrals  each  over  a  cylindrical  volume; 


1  =  f  (1  "  <t)dV  =  2V  Ayi  f6  0  -  {j“)rdr 


(D-6) 


1  =  2tt  11^%  f  o  -  j-)  (f-)d  (f-j 

i=l  Jo  “  e  re 

Now,  break  the  integral  over  r/rg  into  N  equal  zones  (each  jth  voli 
having  the  same  number  of  zones), 
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Finally,  from  (D-5)  and  (D-7), 


('  - » 


M  ,M 

E  j(' 

i=l  j-1 


3  i=l 


The  10.4  diameter  rise  SF6/N2  data  were  used  for  this  example.  A  total 
of  eight  cuts  (at  y  =  2/3  yr  1/3  yT>  central,  1/3  yg,  2/3  yg  and  wake 
cuts  at  3/10  D,  1/2  D,  and  3/4  D)  along  which  - —  were  known  were  used  in 
the  calculation.  The  following  geometry  was  used: 


re2  1— - gr - 

e7l~ - - 


2/3  yj 
V3  yT 

CENTRAL  CUT 


The  values  of  re .  and  Ay..  are  given  below: 


re  = 

f  r(2/3  yT) 

> 

<< 

M 

II 

1/3  yT 

Y'  = 

r ( 2/3  yT)  +  r(l/3  yT) 

Ay2  = 

1/3  yT 

'  e 

2 

Y '  S 

■"(1/3  yT)  +  r(CC) 

Ay3  = 

1/3  yT 

■  P 

e3 

2 

>"4  = 

r(CC)  +  r( 1/3  yg) 

2 

Ay4  = 

1/3  yB 

rP  = 

>"(1/3  yB)  +  r ( 2/ 3  yg) 

Ay5  = 

1/3  yB 

r  5 

2 

r6  = 

■"(2/3  yB)  +  r(0.3  D) 

2 

Ay6  = 

0.3  D  -  2/3  y. 

r7  = 

r(0.3  D)  +  r(l/2  D) 

2 

Ay  7  = 

0.2  D 

*"8  = 

r(3/4  D)  +  r(l/2  D) 

2 

Ay8  = 

.25  D 

also,  r(CC)  =  1/2  D. 

The  density  field  at  the  middle  of  each  zone  is  approximated  by  the  average 
of  the  known  density  field  along  its  upper  and  lower  boundaries  of  each 
zone  (see  Appendix  E).  The  density  in  the  top  zone  is  approximated  as  4/5 
of  the  density  along  the  2/3  y^  cut. 

The  total  mass  was  calculated  to  be  86.4  percent  of  the  initial 
mass  which,  taking  into  account  the  approximations  made  for  this  calculation, 
augments  the  credibility  of  Figure  23.  Furthermore,  the  relative  amounts 
of  mass  within  each  zone  were: 
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Jk 


A  M  11  -'•-a 


Zone 

1 

.0498 

Zone 

2 

.1318 

Zone 

3 

.2009 

Zone 

4 

.1986 

Zone 

5 

.1564 

Zone 

6 

.1695 

Zone 

7 

.0535  j 

Zone 

8 

.0395  | 

Fireball  Wake 


These  results  indicate  that  only  9.3  percent  of  the  total  mass  calculated 
within  the  cloud  is  present  in  the  wake. 


D-2.  MASS  OF  ENTRAINED  AMBIENT  GAS  IN  FIREBALL 

Having  the  mass  of  nitrogen  in  the  mixed  gas  region,  it  was  a  simple 
matter  to  calculate  the  mass  of  entrained  SF6. 

The  mass  of  entrained  SF6  in  the  flow  field  is  just  the  total  mass 
less  the  mass  of  nitrogen; 


mS  *  mT  -  mN  (D-9) 

The  total  mass  is  the  density,  p,  integrated  over  the  volume, 

mT  f  Pdv  =  P.  /  ^  dv  =  P„  |  f  dV  -  f  (1  -  £-)dv|  (D-l 0) 

The  integral  /  (1  -  )dV  was  calculated  in  Appendix  D. 

^  p» 

The  volume  is  simply 


f  M 

J  dV  =  11  £ 

i=l 


rii*y i 


(D-l  1 


Thus,  mT  is  known.  Knowing  mN  and  mN  from  D-l,  the  ratio  of  the  entrainec 
mass  to  the  initial  low  density  mass  may  be  calculated. 

The  result  was  that  approximately  600  times  the  initial  mass  was 
entrained  at  the  time  the  center  of  the  event  rose  to  10.4  initial  bubble 
diameters . 
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APPENDIX  E 


TABULATED  MEAN  CONCENTRATION  DATA 


P  (Atm) 

Gases 

t  (sec) 

D0  (in.) 

2'D0 

D/D0 

Re 

Cut 

Page 

1/2 

Air/ He 

0.48 

1 .0 

5.01 

2.88 

650 

Central 

83 

1 

Air/ He 

0.19  ; 

0.9 

3.16 

2.22 

1100 

Central 

84 

1 

Air/ He 

0.19 

0.9 

3.16 

2.22 

1100 

1/2  D  Wake 

85 

1 

Air/ He 

0.19 

0.9 

3.16 

2.22 

1100 

3/4  D  Wake 

86 

1 

Air/He 

0.19 

0.9 

3.16 

2.22 

1100 

D  Wake 

87 

1 

Air/ He 

0.55 

1 .0 

5.92 

2.85 

1300 

Central 

88 

4 

Air/He 

0.80 

1 .0 

5.65 

3.12 

5200 

Central 

89 

8 

Air/ He 

1 .44 

1 .0 

9.44 

4.84 

10,400 

Central 

90 

10 

Air/ He 

0.90 

1 .0 

5.76 

4.04 

13,000 

Central 

91 

10 

sf6/n2 

0.68 

0.5 

5.99 

3.80 

27,600 

Central 

92 

10 

sf6/n2 

1  .18 

0.5 

1C. 39 

5.12 

27,600 

Central 

93 

10 

SFs/N2 

1.18 

0.5 

10.39 

5.12 

27,600 

2/3  yT 

94 

10 

sf6/n2 

1.18 

0.5 

10.39 

5.12 

27,600 

1/3  yT 

95 

10 

sf6/n2 

1.18 

0.5 

10.39 

5.12 

27,600 

1/3  yB 

96 

10 

sf6/n2 

1  .18 

0.5 

10.39 

5.12 

27,600 

2/3  b 

97 

10 

sf6/n2 

1 .18 

0.5 

10.39 

5.12 

27,600 

0.3  D  Wake 

98 

10 

sf6/n2 

1 .18 

0.5 

10.39 

5.12 

27,600 

1/2  D  Wake 

99 

10 

sf6/n2 

1 .18 

0.5 

10.39 

5.12 

27,600 

3/4  D  Wake 

100 

10 

sf6/n2 

2.50 

0.5 

16.9 

7.26 

27,600 

Central 

101 
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\  r 


P  -  1/2  Atm 
Air/ He 


Re  'v*  650 
Central  Cut 


t  =  0.48  Sec 
DQ  =  1.0  In. 


z/D0  =  5.01 
D/DQ  =2.88 


0.19  Sec 
0.9  In. 


Re  *  1100 
Central  Cut 


z/Dq  =3.16 
D/DQ  =2.22 


P/Pco 

r  (in.) 

r/rebge 

CM 

p/p. 

9.709E-1 

.5 

.5 

4.G07E-2 

7.997E-1 

9.656 

.52 

.52 

3.960 

8.016 

9.589 

.54 

.54 

3.884 

8.046 

9.613 

.56 

.56 

3.844 

8.063 

9.617 

:  .58 

.58 

3.827 

8.070 

9.574 

.6 

.6 

3.798 

8.082 

9.531 

1  .62 

.62 

3.682 

8.129 

9.496 

.64 

.64 

3.674 

8.144 

9.415 

.66 

.66 

3.565 

8.178 

.68 

.68 

3.456 

8.224 

9.276 

.7 

.7 

3.335 

8.275 

9.212 

.72 

.72 

3.191 

8.337 

9.113 

.74 

3.061 

8.394 

EK 

.76 

.76 

2.913 

8.460 

8.331 

.78 

.78 

2.727 

8.544 

8.819 

.8 

.8 

2.620 

8.593 

8.722 

.82 

.82 

2.485 

8.655 

.84 

.84 

2.264 

8.760 

8.405 

.86 

.86 

2.085 

8.847 

8.254 

.88 

.88 

1  .829 

8.974 

KESH 

.9 

.9 

1 .619 

9.081 

7.965 

.92 

.92 

1 .414 

9.188 

7.987 

.94 

.94 

1 .220 

9.292 

7.965 

.96 

.96 

9.73E-3  ! 

9.427 

7.969 

.98 

.98 

6.69 

9.599 

90 


P  =  1  Atm  t  =  0.19  Sec  z/Dq  =  3.16  Re  -v  1100 

Air/He  Dq  =  0.9  In.  D/DQ  =2.22  D  Wake  Cut 


■  . .  - — - 


P  =  1  Atm 
Air/He 


0.55  Sec 
1.0  In. 


D/DQ  =2.85 


Re  1300 
Central  Cut 


r  (in.)  r/r 


r  (in.)  r/r 


CM 

8 

CL 

Q. 

9.54E-3 

9.437E-1 

9.42 

9.444 

9.26 

9.453 

9.10 

9.462 

8.97 

9.469 

8.62 

9.489 

8.35 

9.504 

8.05 

9.521 

7.73 

9.539 

7.45 

9.555 

7.13 

9.573 

6.76 

9.595 

6.35 

9.618 

5.91 

9.644 

5.42 

9.673 

4.98 

9.698 

4.63 

9.719 

4.19 

9.745 

3.65 

9.777 

3.19 

9.804 

2.81 

9.828 

2.55 

9.843 

2.25 

9.861 

1.87 

9.884 

0. 

1  .  OF-^O 

_ _l 

P  =  4  Atm 
Ai r/He 


t  =  0.80  Sec 


Re  *  5200 
Central  Cue 


P  =  8  Atm 
Air/ He 


Re  a,  10,400 
Central  Cut 


T 


t  =  1.44  Sec 
D0  =  1.0  In. 


r  (in.) 

r/redge 

CM 

P/Pm 

0. 

0. 

3.71E-4 

9.977E-1 

4.84E-2 

.02 

4.61 

9.971 

9.68 

.04 

2.21 

9.986 

1  45E-1 

.06 

-1 .83E-5 

1 .000E+0 

1.94 

.08 

6.26 

9.996E-1 

2.42 

.1 

1.01E-4 

9.994 

2.90 

.12 

1.83 

9.989 

3.39 

.14 

1.61 

9.990 

3.88 

.16 

2.75 

9.983 

4.36 

.18 

3.38 

9.979 

4.84 

.2 

2.88 

9.932 

5.32 

.22 

2.40 

9.985 

5.80 

.24 

2.28 

9.986 

6.29 

.26 

3.12 

9.981 

6.78 

.28 

2.76 

9.983 

7.27 

.3 

3.07 

9.981 

7.76 

.32 

3.72 

9.977 

8.24 

.34 

3.61 

9.977 

8.72 

.36 

4.03 

9.975 

9.20 

.38 

5.53 

9.966 

9.63 

.4 

6.47 

9.960 

1 .016E+0 

.42 

8.45 

9.947 

1.064 

.44 

1 .05E-3 

9.935 

1.112 

.46 

1.35 

9.916 

1 .160 

.48 

1.64 

9.898 

z/DQ  =  9.44 
D/DQ  =4.84 


r  (in.) 

r//|  edge 

CM 

P/Pco 

1 .209E+0 

.5 

2.18E-3 

9.865E-1 

1.258 

.52 

2.86 

9.825 

1.307 

.54 

3.41 

9.791 

1.356 

.56 

3.59 

9.780 

1 .405 

.58 

3.61 

9.779 

1.454 

.6 

3.57 

9.782 

1 .503 

.62 

3.55 

9.783 

1 .552 

.64 

3.47 

9.788 

1 .600 

.66 

3.42 

9.791 

1 .648 

.68 

3.26 

9.800 

1 .696 

.7 

3.13 

9.808 

1.744 

.72 

3.01 

9.816 

1.792 

.74 

2.82 

9.827 

1.840 

.76 

2.76 

9.830 

1 .888 

.78 

2.55 

9.843 

1.936 

.8 

2.45 

9.849 

1 .984 

.82 

2.25 

9.861 

2.032 

.8' 

2.21 

9.864 

2.080 

.86 

1.62 

9.900 

2.128 

.88 

1  .71 

9.894 

2.176 

.9 

1 .56 

9.904 

2.224 

.92 

1.36 

9.915 

2.272 

.94 

1.17 

9.927 

2.320 

.96 

9.53E-4 

9.941 

2.369 

.98 

6.47E-4 

9.960 
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0.90  Sec 


Re  'u  1  3,000 
Central  Cut 


4.04E-2 
8.08 
1.21E-1 
1  .62 
2.02 
2.42 
2.83 
3.23 
3.64 
4.04 
4.44 
4.85 
5.25 
5.66 
6.06 
6.46 
6.87 
7.27 
7.68 
8.08 
8.48 
8.89 
9.29 
9.70 
1 .01 0E+0 


.051 E+0 


.091 

.131 

.172 

.212 

.253 

.293 

.333 

.374 


1.4  14 
1.455 
1  .495 
1  .535 
1.576 
1.616 
1.657 
1 .697 
1.737 
1.778 
1.818 
1 .859 
1 .899 
1  .939 
1 .980 
2.020 


5.55E-3 

9. 

5.27 

9. 

4.99 

9. 

4.75 

9. 

4.49 

9 

4.31 

9 

4.09 

9 

3.95 

9 

3.96 

9 

3.84 

9 

3.74 

9 

3.80 

9 

3.69 

9 

3.27 

9 

2.91 

9 

2.84 

9 

2.79 

9 

2.60 

9 

2.36 

9 

2.01 

9 

1  .67 

9 

1.37 

9 

1.15 

9 

8.66E-4 

9 

p  =  10  Atm 

sf6/n2 


.99 

.80 


Re  %  27,600 
Central  Cut 


t  =  0.68  Sec  z/Dq  =  5 

DQ  =  0.5  In.  D/Dq  =  3 


98 


10  Atm 


SF  c/N; 


r  (in.)  r/r 


1 . 24E-1 
1  .48 


t  =  1 .18  Sec 


DQ  =  0.5  In. 


2.05E-3  I  1 .009E+0 


1 .88 
1 . 3GE-4 


1  .007 


1.005 


9.24E-4  1.004 


1 .003 


1.002 


1 .001 


1 .001 


-7.26E-5  1.000 


9.999E-1 


9.995 


9.990 


9.987 


9.980 


9.973 


9.966 


1.02E-3  9.957 


9.946 


9  935 


9.920 


9.893 


9.877 


9.880 


9.882 


9.885 


2/Dq  =  10.39 


Re  ^  27,600 


D/Dq  =  5.12 


6.42E-1 


1 .01 3E+0 
1  .037 


1  .062 


1  .087 


1.112 


1  .136 


1  .161 


1 .186 


1.210 


1.235 


1/3  yT  Cut 


r  (in.)  r/r 


2.72E-3 

2.66 


9.35E-4 


9.887E-1 


9.889 

9.892 


9.897 

9.902 

9.906 

9.908 


9.912 

9.914 


9.915 


9.918 


9.920 


9.921 


9.924 


9.926 

9.928 

9.930 

9.935 


9.941 


9.947 

9.952 


9.957 

9.961 


9.969 


1 .0E+1 


II 


P  =  10  Atm 

sf6/n2 


t  =  1.18  Sec 

=  0.5  In. 
o 


z/DQ  =  10.39 
D/Dq  =  5.12 


Re  *  27,600 
0.3  D  Wake  Cut 


r  (in.) 

r/,redge 

CM 

P/P. 

r  (in.) 

r ' redge 

CM 

M 

8 

Q. 

Q. 

2.21E-2 

.02 

1 .11E-4 

9.995E-1 

5.76E-1 

.52 

9.87E-4 

9.957E-1 

4.43 

.04 

3.56 

9.985 

5.98 

.54 

1.023E-3 

9.957 

6.64 

.06 

4.55 

9.981 

6.20 

.56 

1  .022 

9.956 

8.86 

.08 

5.12 

9.979 

6.42 

.58 

1 .033 

9.957 

1.11E-1 

.1 

5.77 

9.976 

6.65 

.6 

1 .039 

9.956 

1.33 

.12 

5.18 

9.978 

6.87 

.62 

1.031 

9.957 

1.55 

.14 

5.70 

9.976 

7.09 

.64 

9.80E-4 

9.959 

1.77 

.16 

6.16 

9.974 

7.31 

.66 

9.44 

9.960 

1 .99 

.18 

5.80 

9.975 

7.53 

.68 

9.38 

9.961 

2.21 

.2 

5.94 

9.975 

7.75 

.7 

9.21 

9.961 

2.44 

.22 

6.35 

9.973 

7.97 

.72 

9.12 

9.962 

2.66 

.24 

7.08 

9.970 

8.19 

.74 

8.97 

9.962 

2.88 

.26 

7.59 

9.968 

8.42 

.76 

8.74 

9.963 

3.1 

.28 

7.99 

9.966 

8.64 

.78 

8.53 

9.964 

3.32 

.3 

8.29 

9 .965 

8.86 

.8 

8.33 

9.965 

3.54 

.32 

8.31 

9.964 

9.08 

.82 

8.19 

9.966 

3.76 

.34 

8.48 

9.963 

9.30 

.84 

7.87 

9.967 

3.99 

.36 

8.88 

9.961 

9.52 

.86 

7.32 

9.969 

4.21 

.38 

9.18 

9.961 

9.74 

.88 

6.79 

9.971 

4.43 

.4 

9.34 

9.960 

9.97 

.9 

6.11 

9.974 

4.65 

.42 

9.63 

9.959 

1 .019E+0 

.92 

5.30 

9.978 

4.87 

.44 

9.81 

9.960 

1 .041 

.94 

4.49 

9.981 

5  09 

.46 

9.56 

9.960 

1 .063 

.96 

3.63 

9.985 

5.32 

.48 

9.44 

9.959 

1 .085 

.98 

2.55 

9.989 

5.54 

.5 

9.66 

9.957 

1.107 

1 .0 

0. 
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P  =  10  Atm 
SF6/N2 


t  =  1.18  Sec  z/DQ  *  10.39  Re  *  27,600 

0Q  =  0.5  In.  D/DQ  =  5.12  3/4  D  Wake  Cut 


r  (in.) 

r^redge 

CM 

p/pm 

r  (in.) 

r^redge 

CM 

P/Pco 

1.75E-2 

.02 

4.28E-4 

9.982E-1 

4.54E-1 

.52 

6.54E-4 

9.973E-1 

3.49 

.04 

4.15 

9.983 

4.72 

.54 

6.47 

9.973 

5.24 

.06 

1.96 

9.992 

4.89 

.56 

6.38 

9.973 

6.99 

.08 

6.20 

9.974 

5.07 

.58 

6.35 

9.973 

8.74 

.1 

6.31 

9.973 

5.24 

.6 

6,09 

9.974 

1.05E-1 

.12 

6.41 

9.973 

5.42 

.62 

5.96 

9.975 

1.22 

.14 

6.67 

9.972 

5.59 

.64 

5.69 

9.976 

1.40 

.16 

6.75 

9.972 

5.77 

.66 

5.56 

9.977 

1 .57 

.18 

7.16 

9.970 

5.94 

.68 

5.42 

9.977 

1.75 

.2 

7.78 

9.967 

6.12 

.7 

5.29 

9.978 

1.92 

.22 

9.09 

9.962 

6.29 

.72 

5.26 

9.978 

2.10 

.24 

6.61 

9.972 

6.46 

.74 

5.16 

9.978 

2.27 

.26 

6.74 

9.972 

6.64 

.76 

5.03 

9.979 

2.45 

.28 

7.10 

9.970 

6.81 

.78 

4.88 

9.980 

2.62 

.3 

7.49 

9.969 

6.99 

.8 

4.72 

9.980 

2.80 

.32 

8.24 

9.965 

7.16 

.82 

4.47 

9.981 

2.97 

.34 

8.50 

9.964 

7.34 

.84 

4.16 

9.983 

3.14 

.36 

8.67 

9.964 

7.51 

.86 

3.88 

9.984 

3.32 

.38 

8.39 

9.965 

7.69 

.88 

3.60 

9.985 

3.49 

.4 

8.00 

9.966 

7.86 

.9 

3.26 

9.986 

3.67 

.42 

7.70 

9.968 

8.04 

.92 

2.89 

9.988 

3  84 

.44 

7.25 

9.970 

8.21 

.04 

2.48 

9.990 

4.02 

.46 

7.16 

9.970 

8.39 

.96 

2.00 

9.992 

4.19 

.48 

6.91 

9.971 

8.56 

.98 

1.41 

9.994 

4.37 

.5 

6.71 

9.972 

8.74 

1.0 

0. 

1 . 0E+0 
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.. - - -  - - -  ......  j, 


r 


P  =  10  Atm 

sf6/n2 


t  =  2.5  Sec 


D  =  0.5  In. 
o 


z/Dq  =  16.9 


D/DQ  =  7.26 


Re  *  27,600 
Central  Cut 


r  (in.) 

r^redge 

cM 

M 

r  (in.) 

r/redge 

CM 

P/Pco 

3.58E-2 

.02 

-1.48E-3 

1 .006E+0 

9.32E-1 

.52 

1 .66E-3 

9.931E-1 

7.17 

.04 

-1.19 

1 .005 

9.68 

.54 

2.07 

9.914 

1 .08E-1 

.06 

9.94E-4 

1 .004 

1 .004E+0 

.56 

2.08 

9.913 

1 .43 

.08 

-7.45 

1  .003 

1  .040 

.58 

2.03 

9.915 

1.79 

.1 

-5.14 

1.002 

1  .075 

.6 

1 .97 

9.917 

2.15 

.12 

-3.24 

1  .001 

1  .111 

.62 

1.93 

9.919 

2.51 

.14 

-1.92 

1.001 

1 .147 

.64 

1 .88 

9  921 

2.87 

.16 

-5.14E-5 

1  .000 

1 .183 

.66 

i  .83 

9.924 

3.23 

.18 

5.48E-6 

1 .000 

1  .219 

.68 

1.76 

9.926 

3.58 

.2 

1.20E-5 

1 .000 

1.255 

.7 

1.71 

9.929 

3.94 

.22 

5.68 

9.998E-1 

1.291 

.72 

1.62 

9.932 

4.30 

.24 

1.17E-4 

9.995 

1.326 

.74 

1  .46 

9.939 

4.66 

.26 

1  .58 

9.993 

1.362 

.76 

1.28 

9.946 

5.02 

.28 

1.96 

9.992 

1  .398 

.78 

1.14 

9.952 

5.38 

.3 

2.b4 

9.989 

1  .434 

.8 

1.05 

9.956 

5.74 

.32 

3.06 

9.987 

1  .470 

.82 

9.67E-4 

9.959 

6.09 

.34 

3.38 

9.986  ; 

1  .506 

.84 

8.67 

9.964 

6.45 

.36 

3.39 

9.986 

1  .542 

.86 

7.86 

9.967 

6.81 

.38 

3.07 

9.987 

1  .577 

.88 

7.11 

9.970 

7.17 

.4 

3.94 

9.983 

1  .613 

.9 

6.40 

9.973 

7.53 

.42 

5.47 

9.977 

1  .649 

.92 

5.67 

9.976 

7.29 

.44 

6.94 

9.971 

1  .685 

.94 

4.87 

9.980 

8.25 

.46 

8.71 

9.963 

1.721 

.96 

3.94 

9.983 

8.60 

.48 

1 .03E-3 

9.957 

1.757 

.98 

2.76 

9.988 

8.96 

.5 

1.22 

9.949 

1  .792 

1 .0 

0. 

1 .0E+1 
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